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Abstract 

Our calculation of the total cross section for inclusive production of tt pairs in hadron colli- 
sions is presented. The principal ingredient of the calculation is resummation of the universal 
leading-logarithm effects of gluon radiation to all orders in the quantum chromodynamics 
coupling strength, restricted to the region of phase space that is demonstrably perturbative. 
We derive the perturbative regime of the resummed series, starting from the principal- value 

Ph. 

resummation approach, and we isolate the perturbative domain in both moment space and, 

i-Q ■ 

upon inversion of the corresponding Mellin transform, in momentum space. We show that our 
perturbative result does not depend on the manner non-perturbative or infrared effects are 
handled in principal- value resummation. We treat both the quark- ant iquark and gluon-gluon 
production channels consistently in the MS factorization scheme. We compare our method 
and results with other resummation methods that rely on the choice of infrared cutoffs. We 
derive the renormalization/factorization scale dependence of our resummed cross section, 
and we discuss factorization scheme dependence and remaining theoretical uncertainties, in- 
cluding estimates of possible non-perturbative contributions. We include the full content of 
the exact next-to-leading order calculation in obtaining our final results. We present pre- 
dictions of the physical cross section as a function of top quark mass in proton-antiproton 
reactions at center-of-mass energies of 1.8 and 2.0 TeV. We also provide the differential cross 
section as a function of the parton-parton subenergy. 



1 Introduction 



In hadron interactions at collider energies, the main production mechanisms for inclusive 
top-antitop quark (tt) production, as modeled in perturbative quantum chromodynamics 
(pQCD), involve parton-parton collisions. The perturbative series begins at second order, 
0{a 2 s ), in the strong coupling strength, a s . At this order, the parton subprocesses are 

i+j^t + t, (1) 

where the initial partons i,j are either a quark- ant iquark (qq) or a gluon pair (gg). In 
higher orders, gluons are radiated in these two production channels, and there are additional 
production channels, such as qg — > ttq. The gluonic radiative corrections to the lowest- 
order channels create large enhancements of the partonic cross sections near the top pair 
production threshold]!], 0- The magnitude of the 0(a^) corrections implies that fixed-order 
perturbation theory will not necessarily provide reliable quantitative predictions for (tt) pair 
production at Fermilab Tevatron energies. A resummation of the effects of gluon radiation 
to all orders in perturbation theory is called for in order to improve the reliability of the 
theory. This was the main motivation for the published resummation calculations of tt 
production|], |, §. 

In a prior paper, we presented a brief exposition of our method of resummation and its 
application to tt production Our guiding principle is the all-orders resummation of the 
universal leading-logarithm effects of initial-state radiation, restricted to the region of phase 
space which is manifestly perturbative. Our purpose in this paper is to expand upon our 
earlier theoretical and phenomenological discussion. We provide a detailed exposition of our 
approach, highlight the similarities and differences with Ref. |3|, ||, and show that our predic- 
tions are independent of the particular regularization of the infrared and non-perturbative 
regions imposed by principal- value resummation (PVR). We present predictions of the inclu- 
sive top quark production cross section as a function of top quark mass in proton-antiproton 

reactions at center-of-mass energies v^S* = 1.8 and 2.0 TeV, and in proton-proton reactions 
at the energies of the CERN Large Hadron Collider (LHC). We demonstrate that the renor- 
malization/factorization scale dependence of our resummed cross sections is very modest. 
At the end, we speculate on modeling the unknown dynamics of the non-perturbative region 
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and its possible contribution to the cross section. 

In "infrared safe" processes, such as the total cross section for e + e~-annihilation to 
hadrons, there is only one physical scale that characterizes the perturbative cross section, 
namely the total energy of the electroweak initial state. In hadron-hadron scattering, accord- 
ing to the factorization theorem of pQCD, the cross sections for the partonic sub-processes 
are the main object of theoretical calculations. In most hard-scattering processes, however, 
partonic cross sections are not free of singularities. Initial-state hadronic interactions require 
that a non-trivial "mass-factorization" be performed in order to absorb a singular part of 
these cross sections into a redefinition of the long-distance parton distribution functions in 
a universal, process-independent way. Once mass-factorization is performed, the object of 
theoretical interest is the short-distance part of the partonic cross sections (the "hard part"). 
The hard part usually depends on more than one momentum scale, typically because of gluon 
radiation. The hard part is calculated as a series in the strong coupling strength of pQCD, 
and the domain of applicability of the series is a function of all pertinent momentum scales. 
In determining the region of applicability of this series, one encounters difficulties including, 
but not limited to, infrared renormalons||. Most discussions (not necessarily in agreement) 
on the range of validity of pQCD have focused on momentum-scaling properties of Green's 
functions ]7|, |8[], but color factor constants are an equally important aspect of the series and 
should be taken seriously taken into account. Both can contribute to the breakdown of 
pQCD. Although the former have a more "dynamical" appearance because they depend on 
momentum, both create difficulties that are combinatoric in nature. Infrared renormalons 
create factorial growth associated with massless vacuum polarization diagrams, and large 
constants enhance the strength of the partonic interactions according to their corresponding 
fields in color space. Among the many virtues of Ref. || is the demonstration that large 
multiplicative color factors are indeed extremely important in realistic situations. Their 
effects transcend standard operator-product-expansion (OPE) arguments to the effect that 
non-perturbative effects are suppressed by at least A 2 /m 2 or approaches that are usually for- 
mal and model-dependent || . We treat this problem in detail in the present paper from the 
realistic viewpoint of applying resummation to a specific reaction of significant phenomeno- 
logical importance. From an operator point of view, however, this is an issue that deserves 
further investigation, to be addressed in a future publication 0. 
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Top quark production, Eq. ([I]), embodies all the issues mentioned above. An instructive 
illustration can be found in the 0(oP s ) calculation of the partonic hard part, or loosely 
speaking, the "partonic cross section" for heavy flavor production |]J. Figures 12 and 13 of 
the first paper in Ref. |I| show that the 0{a z s ) hard part is much larger than the leading- 
order O{of) hard part in specific regions of partonic subenergy. Analogous results are shown 
here in Fig. 7, to be discussed in section 5. Not all of the range of the partonic subenergy 
can be considered part of the perturbative domain despite the fact that, at finite orders, the 
resulting radiative corrections are integrable throughout that phase space and yield finite 
inclusive predictions. Dominance of the higher-order contributions in the qq and gg channels 
near threshold, and in the gg channel at large values of the subenergy, feeds back (albeit not 
as enhanced) to the physical cross sections obtained after convolution with parton densities. 
In this paper, we limit our attention to processes in which the near-threshold region is the 
most important influence. At Fermilab Tevatron collider energies, the next-to-leading order 
enhancement of the top quark cross section is approximately 25% relative to the leading-order 
value. For comparison, at a scale specified by the mass m of the top quark, a s (m) ~ 10%. 
Top quark production at the Tevatron is one that involves multiple QCD scales. 

Dominance of the next-to-leading order contributions to top quark production near 
threshold in the qq and gg channels is mentioned earlier in this Introduction as a primary 
motivation for theoretical study of resummation of the effects of gluon radiation to all orders. 
One method of resummation was implemented for ti pair production some time ago[0, ||. In 
most resummation methods, the threshold corrections are exponentiated into a function of 
the QCD running coupling strength, a s , evaluated at a variable momentum scale which is 
a measure of the radiated momentum. To recover the inclusive cross section, one then inte- 
grates over the available radiation phase space down to the partonic threshold. An inherent 
ambiguity and limitation of the method of Ref. || is the introduction of undetermined in- 
frared (IR) cutoffs. The need for these cutoffs arises because the kinematically allowed region 
meets the Landau pole of the QCD running coupling. More generally, the exponentiation of 
logarithmic corrections results in essential singularities, whether expressible as singularities 
of the running coupling or not, while a finite-order expansion in terms of a s at a fixed hard 
scale exhibits a polynomial (and hence integrable) dependence on the logarithms. Since the 
IR cutoff dependence is exponentiated in Ref. [0] , the sensitivity of their predicted cross sec- 
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tions to the value of the cutoffs is very significant numerically, especially when multiplicative 
color factors are large, as in the gg production channel. 

An advantage of our resummation approach [[J is that it does not require arbitrary IR 
cutoffs. We use the principal- value resummation (PVR) technique pOj] to by-pass the Landau 
poles and associated renormalon singularities, and we obtain a mathematically unambiguous 
expression for the resummation exponent in moment space. The PVR procedure allows us to 
identify the perturbative properties of the resummation exponent and to separate these from 
the non-perturbative behavior of the exponent. The perturbative properties are obtained 
through a short-distance asymptotic approximation valid in a specific region of moment 
space. The separation is in turn used to derive a perturbative regime for the hard part itself 
in moment space, which is more restricted due to large color factors and exponentiation. The 
inversion of the Mellin transform JTT|] provides the hard part directly in momentum space, 
and we obtain a resummed perturbative expression for the partonic cross section, along with 
the corresponding determination of the perturbative regime in momentum space. This new 
method should be contrasted with the resummation approach of Ref. || where explicit use 
of IR cutoffs makes a perturbative separation impossible. 

It may be argued that adoption of the PVR technique is the adoption of a model and that 
the IR cutoff dependence of Ref. || is merely concealed. We argue that this is not the case 
precisely because we apply our method only in a well defined perturbative region. To the 
extent that PVR is a model, it is a model only for the non-perturbative region that we do not 
include. The main physical issue on which we focus is not the finiteness PVR regularization 
imposes. Rather, it is how this finiteness helps us to probe the asymptotic properties of 
the perturbative cross section. Therein lies part of the usefulness of PVR, even though it 
is possible that the regularization it imposes throughout phase space, including the non- 
perturbative regime, may be physically significant. For example, PVR applied to the Drell- 



Yan process at fixed target energies is in excellent agreement with experiment ||12||. We will 
discard the cross section in the non-perturbative region as model-dependent, notwithstanding 
its finiteness in PVR. 

The perturbative regime should be independent of infrared regulators, including the 
principal-value regulator. To demonstrate this independence, we show in detail in this paper 
how one may obtain the same expressions as ours for the perturbative resummed cross 
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section, independently of any regularization. To do so, we must make two physically plausible 
assertions about the perturbative behavior of the resummation exponent in moment space, 
assertions that are natural ingredients of the PVR exponent. 

Our resummation includes the leading large threshold logarithmic contributions to all 
orders in perturbation theory. These contributions are interpreted to be universal, charac- 
teristic of the initial states that produce the hard scattering. They are common with other 
hard scattering processes, notably massive lepton-pair (11) production, the Drell-Yan pro- 
cess. A complete 0(a 2 s ) calculation is available for the Drell-Yan process, meaning that the 
logarithmic structure is known explicitly in that case to one order in perturbation theory 
higher than for tt production. Final-state effects, which differentiate among hard-scattering 
processes, produce subleading logarithmic structures that are not universal and should not 
be included in a universal resummation approach. It is possible, and even probable, that 
resummation of subleading large logarithms from final-state and interference effects can be 
accomplished, but such a resummation will be process-dependent. For the application we 
have in mind, its numerical effects will be subleading in nature. We consider the resummation 
of final-state logarithms beyond the scope of this work. 

An outline of this paper is follows. In section 2 we present the kinematics of the inclusive 
cross section for top quark production and the form of its resummed version. We include 
a general discussion on counting powers of logarithms for the hard part at finite orders 
and for its resummed version. In section 3 we present the perturbative properties of the 
PVR resummation exponent, E, in moment space, in both the MS and DIS factorization 
schemes, as well as its renormalizat ion-group (RG) invariance properties. We consider a 
formally identical exponent but without regularization, and we examine its factorial growth 
in moment space in order to compare the resulting perturbative asymptotic approximation 
at fixed moments with the one resulting from PVR. The two perturbative approximations are 
in agreement with each other, but the one without the PVR prescription does not directly 
constrain the moment variable since it derives from an unregularized expression that exists 
only formally. We offer additional physically motivated criteria one could use in this latter 
approach to constrain the moment variable, and we show that the resulting constraints are 
similar to those imposed by PVR. Using a heuristic argument, we describe how one may 
estimate the perturbative regime for the exponentiated hard part itself in moment space. 
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In section 4 we discuss the universality of the threshold logarithmic corrections in com- 
parison with massive lepton-pair (11) production, using threshold asymptotics and the results 
of the complete 0(a 2 s ) calculation of heavy quark production |], §]. We examine the detailed 
structure of the hard-scattering function in the gluon radiation phase space (momentum 
space) and the corresponding determination of the perturbative regime in that space. The 
hard-scattering function is proportional to exp(E), which resums the threshold corrections 
directly. The inversion of the Mellin transform produces a hart part in momentum space that 



includes a series of sub leading structures deriving from ^[IT|. It is this latter expression that 
is best suited for the determination of the perturbative regime in momentum space. The two 
determinations of the perturbative regime for the cross section in moment and momentum 
space are in good agreement with each other. We conclude that the resummed pertur- 
bative partonic cross section derived from the PVR approach is independent of the PVR 
regularization and is valid in a well-established perturbative regime. It is derived basically 
from simultaneous minimization of the factorial growth of IR renormalons and color-factor 
combinatorics in the exponentiated hard part. 

In section 5 we discuss the numerical properties and phenomenological behavior of the 
resummed partonic and physical cross sections, for both the qq and gg production channels, 
in the MS and DIS factorization schemes. Predictions are presented in the form of both 
figures and tables of the inclusive cross section for top quark production as a function of 

top mass in proton-antiproton collisions at = 1.8 and 2.0 TeV. We display the factoriza- 
tion/renormalization scale dependence of the physical cross sections and provide extensive 
comparisons with the corresponding quantities of Ref. ||. We address the issue of pertur- 
bative theoretical uncertainties in both our approach and that of Ref. 0, [| . In section 6 we 
present estimates of non-perturbative uncertainties based on physically motivated assump- 
tions about the behavior of the cross section in that regime. We summarize our conclusions 
in section 7. 

2 Production kinematics and resummation 

In this section we begin with expressions for the partonic and physical cross sections in 
finite orders of QCD perturbation theory and the associated kinematics of ti production. 
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Subsequently, we present a resummed expression for the partonic cross section in PVR, 
based on universality of the leading threshold corrections with those in massive lepton-pair 
ill) production. In this paper, we use upper-case S to denote the square of the total energy 
in the hadron-hadron system and lower-case s for the square of the energy in the partonic 
system. 

2.1 Next-to- leading-order cross section and kinematics 

We start with the next-to-leading-order one-particle inclusive partonic differential cross sec- 
tion in the DIS factorization scheme. We use the notation a(/x) = a s (/i)/7r, where fi is the 
common renormalization/factorization hard scale of the problem. A perturbative quantity 
R is expanded as Rfi = 2!- =0 a; J -R , where is the order a J radiative correction to R, 
above the leading order. Unless otherwise specified, a = a(fi = m) where m is the mass of 
the top quark. Following the notation of Ref. || for the subprocess 

-> tfa) +t(p a ) +g(k), (2) 

and defining the partonic invariants 

s = (ki + k 2 ) 2 , h = (k 2 - P2) 2 - m 2 , m = (h - p 2 ) 2 - m 2 , s A = s + t x + «i, (3) 
we express the partonic differential cross section as 

d 2 a (l) 

m 2 s 2 —^-(s,t 1 ,u 1 ) = 2C ij D 1 (s A -,A)ag(s,t 1 ,u 1 )+R. (4) 

In Eq. (f|), o~jj(s,ti,Ui) is the lowest-order Born cross section, expressed in terms of three- 
particle final-state variables, and the remainder R stands for terms that do not contain 
leading logarithmic corrections near threshold. The quantity CV, is the color factor for the 
ij production channel, and 

D 1 (s 4 ,m 2 ] A) = —\n(^)Q(s i -A)+m 2 6(s i )l\n 2 (^) , (5) 

where A is understood in the distribution sense: 

ry 

hiny d(s 4: /m 2 )(f)(s4/m 2 )Di(s4, m 2 ; A) . (6) 
8 



In Eq. (|]), S4 is a measure of the inelasticity of the radiative process, and it is proportional 
to the gluon momentum k. One can show that S4 = 2k -pi and hence, for soft gluons, S4 — ► 0. 
The distribution D\ can be shown to be a "plus" distribution. Indeed, define s^/m 2 = 1 — z 
and 

Then, for any smooth function 0, 

' dzMl - z)D 1 (z) = I dzWl -z)- 0(0)1 ^j 1 "^ - 0(0) rig 1 *' 1 ~ Z > . (8) 

l-j/ Jl-y 1 — Z JO 1 — Z 

The right-hand-side can be written as 

1_£ , ,/ xln(l -2) (/ N r 1 " 6 , ln(l - z) 
dz<t>{l - z)—; --0(0)/ dz- 



1-y 1 — 2 Jl— y 1 — Z 

+ /* Mi - z)» Ki - z) - 0(0) . (9) 

jJTj k\ Jl-e JO 1 — Z 



The series tends to zero as e — > 0, and the result is 



C dzMl - z)D 1 (z) = f 1 6 dz<p(l - z) 1 -^ — ^ + 0(O)-ln 2 e 

Jl-y Jl-y 1 — Z 2 

f ^0(^m 2 )m(4)+ 0(0)1 In 2 e. (10) 
Je S4 \m z J 2 



This is the result obtained also from Eqs. (0) and (|6[) with the identification e = A/m 2 . Hence 
we have proved the identity Di(s^ m 2 ; A) = Di(z), with s^/m? = 1 — z. The identification 
S4/171 2 = 1 — z is suggestive of the similarity of the present reaction, ti production, with 
the Drell-Yan process |T0|, [12] where z = Q 2 / s is the fraction of the squared invariant energy 



carried by the dilepton pair. The threshold is at z — > 1, as it is in our case as well. We use 
z rather than S4 in this paper to stress the similarity between these two reactions. 
We can write the differential cross section, including the Born term, as 



m 2 5 2 ^^-(g, t u m) = h{l -z)+ a2C ij ^ ) )af j (s,t 1 ,u 1 ) . (11) 



If we compare the above expression and the corresponding one for the Drell-Yan process, 
we may verify that the logarithmic structure is identical through next-to-leading order. It 
is important to note here that the exact calculation of the cross section for ti production, 
including the remainder term R in Eq. (|4|), contains sub leading terms, such as (1/(1 — 2:))+, 
as well as constants. These structures are not common to the Drell-Yan reaction, and they 
cannot be resummed as part of initial-state radiation only. As in Ref. |§, we disregard these 
subleading structures in our resummation. We demonstrate the non-universal character of 
these subleading corrections in section 4. 

Following Ref. || we integrate Eq. flTTf ) over the whole partonic phase space. Using the 
appropriate kinematic bounds, we obtain 

4 ] ( ?7 ,m 2 )=/ 1 dz\s(l-z)+a2C ij ( H y z) ) )a*(mz,m 2 ), (12) 

Ji-4(i+ v )+4*A+n [ \ l-z J + ) J 

where 



D _ \/(s/m 2 — 1 + z) 2 — As/m 2 n „ 
a?(v,z,rn 2 ) = ^11 '— dcoB0**{s, h, m) . (13) 

The kinematic transformations 

ti = ~2 m2 ( y s / m ' 2 — 1 + -2 — \J {s/m 2 — 1 + z) 2 — As/m 2 cos 6*^ 

u\ = ~2 m2 (si 171 " 2 — 1 + z + sj (s/m 2 — 1 + z) 2 — As/m 2 cos 

s = Am 2 (l + 7]) . (14) 



are used to obtain Eqs. (|12D and (0). The lower limit of integration in Eq. (0) is derived 
from the kinematics of Eqs. (0) and (^), namely S4 < s — 2m^/s. 

As in the Drell-Yan case, one may eliminate the delta function and plus-distribution in 
Eq. (|12D by integrating by parts. Using the identity 

/ dx(f(x)) + = - f dxf(x) , (15) 

Jz JO 

and the kinematic constraint 



erg ( ri, z = 1 - 4(1 + 77) + Ayjl + T}), m 2 ) = , (16) 

10 



we arrive at the following expressions for the cross sections. In the DIS factorization scheme, 



4 ] ( V ,m 2 ) = f 1 _dz{l + aC ij \n 2 (l-z)}a' ij (r ] ,z,m 2 ) , (17) 
and, in the MS scheme, 

,.\ 

a\] ] {r],m 2 )= _dz fl + a2C ij ln 2 (l - z)\ aUrj, z, m 2 ) . (18) 

yi-4(l+j7)+Vl+»7 L " ' 

o\j{j},z,m 2 ) = ^o%{j},z,m 2 ) . (19) 
The explicit expression for the derivative of the Born cross section in the qq channel is 



a' q q(r),z,m 2 ) = ^C F ir 3 a 2 ^x(z) i^Jx 2 (z) - At + 



2r 



\Jx 2 (z) - At 



(20) 



where 

X(Z) EE 1 - (1 - Z)T , (21) 

r ee m 2 / s = [4(1 + 77)] -1 , and C F = 4/3. In the gg channel, 



3 T 



4r \ 3 /2 

CA-W- ' 



x 2 (^) . 

I t 24r 2 \ , , , / 4r 32r 2 



/ 4r 24r z \ , , , / 4r 32r z \ 1 



37V * 2 w * 2 w w vs 3^w;yr3^ 

where CU = 3 and 



I + a/1 



4r 



= : V X T ■ ( 23 ) 



4t 



2.2 Resummation and power counting 

The large logarithms near threshold, both at finite-orders and in resummation, play a major 
role in our considerations. In this section, we provide a general description that establishes 
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the notion of leading and subleading logarithmic structures and the way we use this termi- 
nology. Generalizing the notation of Eqs. ( |17D and (|T8|) , we write the partonic cross section 
resulting from a finite-order perturbative calculation as 



and c(l,m) are calculable numerical coefficients. In this representation, the classification of 
the logarithmic structures is obvious. Leading logarithms are the monomials proportional 
to c(2m, to), first subleading logarithms the monomials proportional to c(2m — 1, to), etc. 
Since we are concerned with threshold enhancements, z — > 1 or, equivalently, x z — ► +00, the 
nomenclature is related directly to the numerical importance of the corresponding logarith- 
mic structures in the cross section. 

At finite orders, the logarithmic structures are integrable near threshold, and their con- 
tributions to the cross section are finite, provided the Born cross sections are also integrable 
(as is indeed the case). For purposes of demonstration, let us assume that the Born cross 
section behaves as (1 — z) v , —1 < v. Using 




(24) 



Here 




(25) 




(26) 



and denoting L(rj) = 1 — 4(1 + if) + 4a/1 + rj, we find 





JL(n) 



x\ = a m c(/,m)(-l) 



' lim 



= a m c(/,m)(-l) / lim 




H ( V )K«> m \ V ,a) 



(27) 



where 
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is the "Born" cross section, and 



K^\ V , a) = a m c(l, m)—^ + v) j ln'(l - L( V )) (29) 

is the "K-factor" due to the c(l,m) radiative correction in Eq. (pop. Near threshold, 

lim L(rj) ~ 1 - 2r] . (30) 

While both Ho{rj) and H^' m \rj,a) vanish as rj — > 0, the function K^ l,m \ri,a) diverges as 
the Z-th power of the logarithm. Relative to the lowest order term, its contribution can be 
arbitrarily large in 77 space for a sufficiently high power I, despite the perturbative suppression 
a m . The logarithmic terms will have a significant effect on the physical cross section, obtained 
as the integral over 77 of Eq. ( |27|) convoluted with the parton distributions, especially if the 
hadronic center-of-mass energy is such that the support of this integration emphasizes the 
threshold region. Under such circumstances we may view 77 as a second important physical 
parameter, in addition to m, that is probed phenomenologically. The two-scale nature of 
the problem is evident in that we have one hard scale m and a ratio rj = (s/4m 2 ) — 1 
whose important physical domain is near 0. The perturbative QCD series is required to have 
reliable behavior in both variables, m and rj. 

The objective of resummation is to derive formulas, based on the properties of QCD as 
a field theory, that provide a summation of the various classes of monomials in Eq. (f25|). 
In typical resummation methods, the partial sums of the hard part, Eq. (p5|), are replaced 
with a resummed hard function that contains the numerically important pieces of the partial 
sums, to all orders in pQCD. Most resummations of threshold effects result in exponentiation 
of the large logarithmic contributions. A generic resummation may result in an expression 
such as 

H(z, a) ~ e E{x *' a) . (31) 

The principal content of resummation resides in the exponent E(x z , a), which is typically a 
function of the QCD running coupling strength, integrated through intermediate momentum 
scales. The particular form of this function is fairly process-independent and follows from 



general field-theoretical argument s||13||. For the present, we concentrate on power-counting, 
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ignoring complications specific to particular resummation methods. A large part of the 
remainder of the paper is devoted to the complications. 

Under specific conditions that we analyze in detail later, the resummation exponent can 
be cast into a perturbative form similar to that of Eq. fl25|) , 

N m+1 

£(x z ,a)=5> m ]Te(Z,m):4. (32) 

m=l 1=0 

The precise power structure above follows from first principles. For a resummation method 
to be realistic in practice, the exponent E(x z ,a) should be calculable in a finite number 
of steps. One consequence is that the exponent is not calculable with arbitrary precision, 
implying, in turn, a level of uncertainty in some of the coefficients e(l, m), for specific ranges 
of I, m. The only consistent treatment of this limitation is to calculate the exponent in enough 
detail that the uncalculable coefficients in Eq. (|32|) accompany logarithmic structures x l z that 
are numerically insignificant in the range where the representation Eq. ([32] ) makes sense in 
the first place. As demonstrated explicitly in the Drell-Yan process|Tl[], and as will become 
apparent below, use of the two-loop QCD running coupling strength is necessary in the 
calculation. Roughly speaking, the representation Eq. ( |3"2"D is valid when ax z < 1. This 
inequality determines the limit of calculational accuracy of the exponent. The upshot is that 
the coefficients e(l,m), I G {0,m — 1}, m > 3 are undetermined. 

Because the exponent is the central object in resummation, most of the approximations 
are effected on that quantity. By way of definition we say that we resum leading logarithms 
when all {e(m + l,m)} are determined, and non-leading logarithms when all {e(m, m)} 
are determined in addition. Resummation of leading and non-leading logarithms has been 
done in detail for the Drell-Yan process |12||. It is important to note, however, that our 



definition is not identical to the obvious definition that exists at finite orders. The latter is 
contained in the former. In the next few paragraphs, we establish the relationship between 

the power-counting of logarithms in resummation and the obvious finite-order definition at 

the beginning of this section. 

Leading-logarithm resummation includes all leading logarithms at finite orders. Indeed, 
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m=l 



suppose we have performed leading-logarithm resummation. Then 

N 

E(x z , a) ~ £ a m e(m + 1, m)x™ +1 , (33) 
and the kernel of the hard part is 

H(z,a)= Y[ e ^e {m+ i,m )x T^ (34) 



N 



m=l 



To compare with finite-order pQCD, we must make a Taylor expansion of Eq. ( p4[) in a. 
It is clear that all leading logarithm coefficients in Eq. (pq) , c(2m, m), are obtained from the 
set {e(m + 1, m)} of Eq. (|34l) and, more specifically, from e(2, 1) itself. Consider for example 
the Taylor expansion of the two first terms of the product of Eq. (|34|) : 



^ q k ^x 2 z kl e kl (2,l) ~ a 2fc2 xf 2 e fc2 (3,2) 

h w x ^ w ' (35) 



The general monomial of this product is 

a fc 1+2 fc 2;E 2fc 1 +3fe 2e fc 1 (2 ? i)e fc 2(3 ; 2) 



h\k 2 \ 



(36) 



and the only terms that fit the set described by c(2m, m) are h\ G N, = 0. Therefore, 
c(2m,m) = e m (2, l)/m\, and leading-logarithm resummation includes all finite-order leading 
logarithms, along with a specification of the values of these terms beyond the order in pertur- 
bation theory at which they may have been computed explicitly. It should be remarked that 
it includes more. The product of Eq. (0) generates upon expansion sub leading logarithms 
in finite-order pQCD. For example, k\ G N, k 2 = 1 in Eq. ( |3~6D are of the first-subleading kind 
c(2m — l,m) of Eq. (|25|) . These are not a closed set, i.e., they account only partially for all 
the subleading logarithms resulting from a finite-order calculation. For example, to account 
for all {c(2m — l,m)}, a non-leading exponent in Eq. ( |3"2"D would have to be calculated, i.e., 
all coefficients {e(m, m)} in Eq. ([32]) would in addition have to be specified. Upon expansion, 
such a resummation of finite-order first subleading logarithms would produce a series 

oo 2m 

H{z, a )=Y, am Y, h (h m Vz , (37) 

771=0 Z=7Tt 
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where each coefficient h(l,m) is a product of {e(7,m)}'s. 

We remark that Eq. fl57|) does not determine uniquely all finite-order coefficients {c(l > 
m, m)}, for all values of m, but it leaves a subset undetermined. To illustrate this point, 
consider the leading uncertainty in the resummation exponent. Since the one-loop running 
coupling constant produces all {e(m + l,m)} and part of {e(m, m)}, m e N, and the two- 
loop running coupling constant produces the remaining pieces of {e(m > 2,m)}, the leading 
uncertainty in the exponent is 8o?x\ in realistic resummations. This is the lowest-order in 
a, highest-order in x z undetermined monomial that comes from expanding the three-loop 
running coupling strength in the exponent. 

A product of the form 

£ — w — "S, - w~~ (38) 

would then be present in the resummed hard kernel. The general monomial is 

a k 1 +3k 2x 2(k 1 +k 2 ) e k 1 ^ 1)(5 fe 2 



(39) 



We wish to investigate the degree to which the second series in Eq. ([38|) maximally changes 
the coefficients {h(l, m)} of Eq. (P7|). Setting 2(k\ + kz) = 2(ki + 3^) — n, we see that 
we obtain k 2 = n/4. Hence the minimum n = 4, and the minimum k 2 = 1. Maximum 
uncertainty arises in the coefficients 

h(2m-4,m). (40) 

The coefficients {h(l > m,m)} are affected for high enough m. For example at minimum 

power of a, m = 4, the "diagonal" coefficient /i(4, 4) is affected, at m = 5 the coefficients 

h(6, 5) , h(5, 5) are affected, etc. However, the leading, first-, second- and third-subleading 

logarithms, as defined in finite-order power-counting, are unaffected for any orders. These 

finite-order logarithmic structures are then resummed to all orders, by a resummation of 

leading and non-leading logarithms. It is unlikely that the affected structures are significant 

numerically in such a resummation, once we restrict ourselves to the perturbative regime in 
x z . 
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We complete this section by returning to the reaction of interest, tt production. At 
this point, we invoke universality with the Drell-Yan case. Because the finite-order leading 
logarithms are identical in the tt and 11 cases, we can resum them in it-production with the 
same function we use in the Drell-Yan case. According to Ref. (Tl|, the structure of the 
kernel of the resummed hard part in the Drell-Yan case is 

(E(x z ,a) oo \ 
- r — ^Q^ar,, a) J . (41) 

For simplicity, we drop the channel indices ij; x z = ln(l/(l — z)). We defer to sections 3 
and 4 an explanation of the various functions in this equation. One can therefore write the 
resummed partonic cross section as 



i 

l-4(l+r))+VT+^ 



a i:j (ri,m 2 ) = dzI(z,a)afAr], z,m 2 ) . (42) 

/ 1 A ( 1 I ^ I A /1 I m J 



After integration to get rid of the delta function and plus-distributions, we find 

a ij (r],m 2 ) = [ dzH(z, a)cr',(?7, z,m 2 ) . (43) 

Jl-A(l+ri)+AJl+r, J 



'l-4(l+J7)+4vTH? 

The kernel of the hard part is 

H(z, a) = l+ dxe E ^ a) ]T Qj(x, a) . (44) 

The kernel in Eq. Q44]) depends solely on the resummation exponent E(x,a), either 
explicitly, or through the functions Qj which depend exclusively on E. This exponent, 
in turn, depends on the factorization scheme, i.e., for the process under consideration, on 
Eqs. ([H]) and flT8|). It is to this exponent that we turn our attention in the next section. 

3 The Resummation Exponent 

Invoking universality with the Drell-Yan case, we can express the exponent for tt production 
using the results of Refs. [T3| . By the same token, because of the restrictions on that univer- 
sality we may keep only the pieces of the exponent that are universal in the two cases, i.e., 
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the pieces that reproduce upon expansion the leading logarithmic structures attributable to 
initial-state radiation, exemplified in Eq. (|17D. The details of the inversion of the Mellin 
transform, expressed in Eq. ([44]), may be found in Refs. JIT , |12]| and are analysed further in 
section 4. 



In this section, we present the exponent E(x, a) of Eq. ( fHj) in moment space, where the 
moment n = exp(x). For the Drell-Yan process, the exponent in moment space in the PVR 
approach may be written in either the DIS or the MS factorization scheme |L0j, fT^) . In the 
DIS scheme, 

E(x,a) = - dC-——\ ^ gi [ a (\ m )}+g 2 [ a ((l-C)m)) , (45) 



and, in the MS case, 



E(x,a) = - / d(±— - / gi [ a (Xm)]. (46) 

JP 1 — C, "'(l-C) 2 A 



Here P is a principal- value contour [10]; gi and #2 are functions of the QCD running coupling 
strength a (Am). It is important to note that Eqs. ( pES]) and ([K]) include in general all large 
logarithmic structures in the Drell-Yan case. This exponent is renormalization-group invari- 
ant, by construction. The logarithmic structures generated by this exponent are recovered 
upon expanding the functions gt as perturbative series with respect to the running coupling 
strength and re-expanding the running coupling strengths in terms of the hard-scale cou- 
pling strength a(m)[|T0j. As we argued in the general discussion of section 2.2, and as shown 
explicitly in ]IpJ, a series of truncations can include all large threshold corrections. 

The only necessary ingredients are the two-loop running coupling strength and a few 
terms in the expansions of the gj's. In general, gi generates all the leading threshold cor- 
rections and some non-leading ones, and g<i completes the resummation of numerically im- 
portant non-leading ones. Since we base our resummation on universality, which is valid for 
leading logarithms only, we disregard #2 for most of the rest of this paper, and we use the 
notation g = g 1 . We apply the consistency requirement of resumming leading logarithms 
only by performing the appropriate truncations and neglecting subleading structures when- 
ever they appear. We describe these structures quantitatively in this section. We discuss the 
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renormalization-group (RG) properties of resummation in section 3.4. The process of resum- 
ming leading logarithms only, and the associated truncations, result in an approximation of 
the renormalization-group-invariant exponent of Eq. (|4~5"D by an exponent that is approxi- 
mately renormalization-group invariant. By varying the factorization/renormalization scale 
H within a logical range about its central value \l — m, we obtain a variation of the resulting 
resummed cross section that serves as a quantitative measure of the effects of logarithmic 
structures that are not resummable in the ti process. Variation with \i is the bulk of the 
theoretical uncertainty. 

3.1 Truncation and perturbative representation of the exponent 

Using the truncations 



OO 



g[a(Xm)} 



E 



(47) 



along with 





we write the PVR exponent in the MS scheme as 




(49) 



This integral can be evaluated exactly |]TD|. The result is 




(50) 



where 




(51) 



and the exponential integral is defined by the principal value 




(52) 



19 



Equation ( |49| ) or, equivalently, Eq. ([5(]) has a perturbative asymptotic representation 



E(x,a) ~ E(x,a,N(t)) = g 



N(t)+1 p+1 

1 E « P E^'- 

p=l j=0 



(53) 



This representation is valid in the moment-space interval 



1 < x = Inn < t. 



(54) 



The coefficients are 



-b^\-l)^2"c p+1 ^{p-l)\/j\ 



(55) 



and the constants Ck are obtained from the expansion T(l + z) — Y^T=o c k zh , where T is the 
Euler gamma function. 

The number of perturbative terms N(t) in Eq. (|53|) is obtained by optimizing the asymp- 
totic approximation 



Using Eqs. 



E(x, a) — E(x, a, N(t)) = minimum. 



and (|55|) , we rewrite Eq. fl5"6|) as 



(56) 



E 

m=l 



rrumf 



1 - n) N{t) (-lYo* p+2 C-lV _1 c • 

- . l m £(mt) - J2 E - — ^ 

p=0 



i=o 



J! 



minimum 



(57) 



Equation Q5"7| ) denotes an approximation between two functions of the moment n, and it 
determines the optimum number of perturbative terms N as a function of the parameter 
t. Indeed, as shown numerically below, within the interval of Eq. ([54]), the optimization of 



Eq. (|57|) has a solution that depends on t only, N = N(t). It also can be shown |L2| that in 
the complementary interval 

t < x = Inn < oo, (58) 

the approximation Eq. (^) breaks down for any integer N that is a function of t only. For 
very large n within the interval of Eq. (ESI) the asymptotic approximation 



X 



(59) 
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holds, and it is clearly a non-perturbative one. 

Throughout this paper, we use the two-loop formula for the fixed coupling strength 



a(m) = ^ = - WlnK/A 2 )) 

{ } ~ 7T & 2 ln(m 2 /A 2 ) b\ ln 2 (m 2 /A 2 ) ' K ] 

with 

b 2 = {11C A - 2n/)/12, 6 3 = (34C1 - (10CU + 6C F ) % )/48, (61) 

and number of flavors n f = 5. We set A = 0.158 GeV (the CTEQ3M valued). 

In Fig. 1 we illustrate the validity of the asymptotic approximation for a value of t 
corresponding to m — 175 GeV. In Fig. la we show how N(t) is determined from Eq. (|57D 
for a fixed t and selected parametric values of n. The plot shows that optimization works 
perfectly, and it demonstrates the typical breakdown of the asymptotic approximation as 
N increases beyond N(t). This rise is the exponential rise of the infrared renormalons, the 
p\ growth in the second term of Eq. (]57|) . As long as n is in the interval of Eq. (|54[) , all 



the members of the family in n are optimized at the same N(t), showing that the optimum 
number of perturbative terms is a function of t only. In Fig. lb we plot the function N(t) 
for a range of t relevant to tt production. An excellent numerical approximation is provided 
by the fit 

JV(t)~[t-3/2], (62) 

where the integer part is defined as the closest integer from either direction. (It is amusing 
that this fit suggests a lower limit for the hard scale m. A perturbative series is an improve- 
ment in accuracy if N(t) > 1 which, from the above fit, implies t — 3/2 > 0.5. Using for 
simplicity a one-loop fit of a(m) and A, we deduce that the hard scale m/A > e 2 ~ 8, well 
within our expectations.) 



Equations (|5|) and ( |57| ) suggest a perturbative behavior for the exponent in moment 
space independently of the color factors which reside in In particular, the range of 

validity of the perturbative expression for the exponent, x/t < 1, is obtained by direct 
comparison with the exact principal-value definition, as shown in Fig. 2. (Equations ( |58| ) 
and (|59D are also suggestive, but they constitute an asymptotic limit, valid for x larger than 
t.) This range of validity has the consequence that terms in the exponent of the form a k ln fc n 
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are of order unity, and terms with fewer powers of logarithms, a k ln fc m n, are negligible. This 
explains why resummation is completed in a finite number of steps in the Drell-Yan process, 
as discussed earlier. The same is true here, but, in addition, we discard monomials a k \n h n 
in the exponent because of the restricted universality between the tt and 11 processes. 
The exponent we use in the rest of the paper is the truncation 

N(t)+1 

E(x,a,N) = g {1) £ a p s p x p+ \ (63) 

with the coefficients 

s p = s^i jP = 6T 1 27p(p + 1) • (64) 
We note in passing that the leading logarithm truncation of the exponent, Eqs. fl63|) and 
(|64|), forms a convergent series, 

E(x, a, N) ~ ^x 1 1 + (— ^ ~ X ) ln ( X ~ 2ah2X "> } ' ( 65 ) 



as long as x < l/(2«6 2 ) = t. This, of course, is not true for the exponent of Eq. 
that contains the full subleading logarithmic structures in the Drell-Yan process. We use the 
convergent version of the exponent Eq. fl65|) in special cases to arrive at simplified expressions 
for the perturbative regime in momentum space. 

There are equivalent expressions appropriate for the DIS factorization scheme. Equa- 
tion (f49|) becomes 

C™- 1 - 1 /-i-f d\ 1 



p /-n—l _ i p 

E(x,a)~-agW / d( \ / 
Jp 1 — C h 



1 — C •Ai-C) 2 A l + a6 2 lnA 
o« f 1 ~ (l_ n ) m 1 » (i_ n ) ] 

The perturbative approximation, subject to the truncation requirements discussed previ- 
ously, is 

N(t)+1 N(2t)+1 

E{x,a,N) = gW £ ^s p x" +1 - g^ £ (-) v^ +1 , (67) 



P =i P =i 
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or, in convergent form, 



E(x, a, AO ~ ^x { ( — | l) ln(l - 2ab 2 x) - ( ~4 l) ln(l - ab 2 x)\ . (68) 

b 2 l\2ab 2 x J \ab 2 x J J 



The convergent exponents, Eqs. flog) and (p8|), suggest the perturbative interval of Eq. 
fl54l), independently of PVR. On the other hand, beyond the end point x = l/{2ab 2 ) these 
convergent expressions are ill-defined because the leading term has a branch-point singularity. 

3.2 Regularization-independence of the perturbative exponent 



It is valuable to stress that we can derive the perturbative expressions, Eqs. (|53|) , 
and (|62|), without the PVR prescription, although with less certitude. The analysis of this 
subsection is presented in order to show that our final perturbative results do not depend 
on the specific manner that the infra-red region is regularized in the PVR approach. 

We begin with the unregularized form of Eq. (|49D , i.e., with the integral over x done on 
the real axis: 

E (x,a) = -ag {1) / dx— / — — — . 69 

v ' Jo 1 - x 7(i-x) 2 X l + ab 2 \nX v 7 

We expand the inner integrand as a Taylor series around a. As we will see, part the problem 
now transforms into ignorance of the asymptotic properties of this expansion. Writing 

(i-s)2 A 1 + a6 2 hi A h pt^ 1 



we deduce 

±1 P^ 1 y 



In" i/. (71) 



The upper limit of the summation in Eq. Q7Ip is left undetermined, because, as is made 
evident below, the series represented is only formal, i.e., not convergent. Lack of convergence 
is associated with the Landau pole exhibited by the original integral, Eq. (|69| ). Using the 
identity 

■"^ S3 (£)'»" • (72) 
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and the Stirling approximation for the beta function B(n,rj), 



( a \p ( i \ p+i /-I y 

S3 (as J [B^*) --)=?' g V" c >«- to ' " ■ < 73 > 

we find 

P+i 

.EbO^; a) = g^ 1 ' ^2 Qp s i.p xJ ' ' ^ = ban, (74) 
p=i 3 -=o 

with the coefficients of Eq. fl55|). This expression for ^(x, a) is the same as Eq. fl53"|), the 
only and major difference being that we do not know the asymptotic properties of this series 
in the full range of moments n. The added information is precisely what is furnished by 
PVR, as we saw earlier. Both the function N(t) and the range of validity of the perturbative 
expression, 1 < x = Inn < t, are provided by the principal- value prescription. 

To be more explicit, we examine Eq. ( 74) in some detail. Because the coefficients 



p 



grow factorially, the series does not represent a convergent infinite series for a fixed t and 
n. The factorial growth is precisely the infrared renormalon induced by the existence of the 
Landau-pole singularity of the original integral. Without this factorial, as in the truncated 
expression Eq. (|63|), the resulting infinite series is convergent for a fixed n, but it diverges 
nevertheless at the threshold n — ► oo due to the powers of the logarithm. Since E is the 
exponent of the cross section, the resulting singularities in the cross section would be essential 



singularities. In fact, in Eqs. (|53|) and (pq) there is a trade-off between factorial growth and 
powers of logarithms: greater factorial growth is accompanied by fewer powers of moment 
logarithms (and hence of momentum logarithms in the cross section). 

We conclude that the nature of the series in Eq. ( [74] ) is asymptotic, and we rewrite 
Eq. ([T§ as 

N p+1 

E (x, a) ~ E(x, a, N) = g (1) J2 aP Yl s j,p xj > x = lnn ■ ( 75 ) 

p=l j=0 

Because the original integral is unregulated, the properties and range of validity of this 
asymptotic series are not obvious. In PVR, regularization is incorporated, and, since there 
are no undetermined extra scales, such as the introduction of IR cutoffs, the asymptotic 
properties are determinable fully, as shown in section 3.1 
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For a fixed t and n, one may use the monotonicity behavior of the corresponding partial 
sums to try to determine an upper limit for the number of terms in Eq. (|75[) . This procedure 
is illustrated in Fig. 3a for m = 175 GeV. We note that beyond a certain range of N, the 
exponent increases factorially, a demonstration of both the asymptotic nature of the series 
and of the effect of the IR renormalons. A range of optimum N can be determined where 
the growth of the sum reaches a plateau, before the factorial growth sets in at large N. The 
exponent is fairly flat in this region so the indeterminacy of the optimum N is insignificant 
numerically The plateau is centered around N opt ~ 6 to 7 for a wide range of moments, 
in agreement with the results of the previous subsection. One can make these statements 
quantitative by defining the slope 

S(x, a, N) = dE ^ d ^ — = E(x, a,N + l)- E(x, a, N) . (76) 
For a range of n where there is a plateau, the optimum N can be determined by the equation 



N opt ■ S(x,a,N) 



= minimum . (77) 

x,a= fixed 



In Fig. 3b, we observe that the optimization of the perturbative exponent with this method 
gives the same results as the principal- value method, i.e., 

No*- l~N(t) = [t- 3/2]. (78) 



A second issue is whether one could determine the same range of validity, Eq. ([54]), of 
the representation Eq. ( [75] ) without PVR. This is possible if we impose the supplementary 
requirement that this representation results in a value of N opt that depends only on the 
hard scale m, and not on the moment n. In Fig. 3b we show Nopt, the minimum of the 
slope S(n, N) for various values of n at m = 175 GeV. We observe that N opt is constant in 
approximately the same range of n as in Eq. (0). When n is increased beyond that range, 
the corresponding N op t becomes a function of n. As n is increased, the plateau created by 
the corresponding values of the exponent starts shrinking, and the renormalons move to the 
left of the figure. For values of n close to the boundary of Eq. ([54]), the slope in Eq. ([76]) 
becomes substantial. 
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The two requirements, a plateau in E(x, a, N) and a value of N opt that depends only on 
the hard scale, allow one to make statements concerning the perturbative representation of 
the exponent that are similar to those the PVR approach, albeit with less certitude. The 
perturbative representation of the exponent, Eqs. fl5^), fl54]) , ([55]), and flS^) , is obtained in 
straightforward fashion in PVR, but the analysis of this subsection shows that the result is 
approximately independent of the regularization PVR imposes. The result can be recovered 
by a study of the factorial growth of the partial sums and the reasonable physical requirement 
that the plateau of stability of the exponent depend on the hard scale of the problem only. 

3.3 Perturbative regime of the hard part in moment space 

We have established a perturbative representation for the exponent in moment space, with 
or without principal-value regularization. The conclusions are similar. The exponent has a 
perturbative representation in the range 

1 < x = In ra <t = (79) 

independent of the constant i.e., of the channel-dependent color factors. This is a 
maximum range allowed for the perturbative cross section in moment space. One should 
distinguish, however, between the range in moment space where the exponent has a per- 
turbative representation, and the range where the cross section itself is perturbative. The 
distinction arises because the cross section, which is proportional in moment space to the 
exponential of E(ln n, a, N), is much more sensitive to variations in n. The cross section also 
depends exponentially on the color factors g^ l \ which can be much larger than unity, while 
the perturbative representation of the exponent in the regime of Eq. (|79|) is independent of 
the color factors. 

We address the question of the perturbative regime in momentum space in section 4, 
where we present the inverse Mellin transform that provides the cross section. However, 
one can furnish a heuristic argument in moment space as well. Again, this argument is 
independent of the regularization prescription of the original integral representation of the 
exponent, once we work with its asymptotic perturbative approximation, Eqs. (Q) and (|63|). 

The idea is to regard the plateau of Fig. 3a as a region of perturbative stability in both of 
the variables it depends on, namely N and n. Denote the kernel of the hard part in moment 
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space by 

I(n,a,N) =e E(lnn ' a > N) . (80) 

Then, this kernel is perturbatively stable in an interval of n such that variations around 
the center of the plateau N(t), which does not depend on n in the perturbative regime, 

provide 0(a s ) variations for I(n,a, N), i.e., numerically negligible contributions that are 
not enhanced by threshold effects. This statement is made quantitative by the requirement 

8 N I(n,a,N(t)) = e E{lnn ' a > N{t)) 5 N E(\nn, a, N(t)) = e E{lnn ' a ' N{t)) S(\nn, a, N(t)) < a s , (81) 



where the slope S is defined in Eq. (|76|) . Another way of reaching the condition of Eq. (|3lD 
begins with the order-by-order perturbative hard part written as 

m{n, a) = £ a m s m {n) + ]T a m r m {n) . (82) 

m=0 m=l 

The first sum in Eq. (|82f) contains singular functions of n that produce large threshold 
enhancements and must be resummed, while the second sum contains regular functions 
of n, whose numerical effects are of 0(a s ). The first sum is replaced by our resummed 

I(n, a, N(t)), and the condition in Eq. ( PT| ) is equivalent to a plateau of perturbative stability 
such that variations of the resummed terms around the plateau belong to the second, regular, 
sum of Eq. (||): 

e E(lnn,a,N(t))+8 N E(lnn,a,N(t)) ^ e E(lnn,a,N(t)) + e E[lnn,a,N(t)) ^ 
I 

~ e BQ*w,m) + a m r' m (n), (83) 

m=l 

so that 8jfl{n,a, N(t)) ~ 0(a s ). This is the same criterion as Eq. (|5T]). 

The regular sums in Eqs. (B2T) and (B3|), have been calculated only to first order in ti 



production. This is the reason we consider them at most equal to a s , rather than assum- 
ing definite numerical coefficients, large or small, unsupported by exact calculations. This 
statement is corroborated by the exact one-loop calculation of the tt cross section, and by 
the one- and two-loop calculations of the Drell-Yan cross section. It is in general assumed 
in all resummation procedures. 
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3.4 Renormalization-group properties of the perturbative expo- 
nent 

We have studied the exponent under the simplification that the renormalization/factorization 
scale, the "hard scale" , is fixed at p = m. In this subsection, we discuss the dependence on 
/i of the resummation exponent in moment space. We work entirely in the MS scheme. Since 
RG invariance is one of the ingredients of resummation, we expect that the full Drell-Yan 
exponent, containing the full scale dependence, is exactly RG invariant. It can be seen by 
inspection that Eq. (BH) is scale- invariant: 



E(x,a) = E(x,a{p),m/p) = - d(— — / — #i[oi(A//)]. (84) 

JP 1 — C, J(l-Q 2 m 2 /n 2 A 

Here we ask to what degree the truncations we impose, starting from Eq. ([49|) , affect 
this invariance. We begin with the explicitly /^-dependent equivalent of Eq. ([49]) , read from 
Eq. fl34|). The corresponding expression isfT4[ 



E(x,a(u),m/u) ~ -a(p)g {l) / dQ- -± / ^ / NI - - , (85) 

where the dependence on /i is both implicit (in ct(/i)) and explicit. Working along the lines 
of section 3.2, we may show that 

E(x, a(/i), m/ p) ~ E(x, a(p),m/ p, N(p)) 

= -gV> £ <» W i)Zb2 EAi'^K^E^W- (86) 
p=i P k=i\P~ K >- j=o 3- 

We use the approximation N(p) ~ N(t(p)). This equation makes clear that explicit p- 
variation (i.e. k < p) is equivalent to inclusion of non-leading logarithms x^,j < p + 1. In 
this sense, uncertainty expressed through /i-variation overlaps with uncertainty expressed 
through inclusion of non- universal logarithms. In the triple sum of Eq. (^Bj), the k = p 
term is the dominant term numerically in the exponent. This term represents the implicit 
/^-dependence, Eq. (0). All other ln(m//i) terms are non-leading in x = Inn. The linear 
term in ln(m//z), obtained for k — p — 1, has a universal Inn substructure. It is the one 
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containing equal powers of a(fi) and Inn. It comes from the same diagrams as the leading 
Inn structures at fx = m and amounts to exponentiation of the explicit 0{oP s ) /^-variation. 
Our /i-dependent exponent for tt production can therefore be written 



£(x,a(u),m/u,AT(u)) = </ (1) | £ a p (fi) ^ +i) xP +1 -\n(m/fi) £ ""M" 



P ) 

(87) 

We call attention to the sign structure of Eq. fl8"T|). For a fixed moment and scales fi decreasing 
from m, the implicit growth of the first sum due to the increase of a(ft) is cancelled by the 
minus sign of the explicit dependence of the second sum, and vice-versa. On the other hand, 
in the truncation from Eq. fl86| ) to Eq. (87), non-universal subleading logarithmic structures 



are discarded. Therefore we expect only approximate scale invariance of Eq. (|87j) . 

In Fig. 4 we plot the normalized exponent, in each of the truncated versions of Eqs. (|S5|), 
fl36|), and fl3?D, versus a wide variation of fi G {100,300} GeV, and for various perturbative 
values of the moment n. We see that both the PV Drell-Yan exponent Eq. fl85|) and its per- 
turbative Drell-Yan version Eq. ( [SB] ) are straight lines, practically coincident. The universal 
exponent Eq. fl5T| ) is a gently increasing function of fx, and it lies somewhat below the Drell- 
Yan exponent. The increase creates a partial compensation of the explicit //-dependence of 
the Born cross sections, more so than if the full Drell-Yan exponent were used. This differ- 
ence stands to reason, since no implicit /i-dependence exists in the Drell-Yan process at the 
Born level. These //-dependence properties conform to the intuition that our resummed cross 
section should show less variation with fi than the next-to- leading order counterpart. We call 
attention to Fig. 4b, where we show the effects of /i-dependence as a function of the number 
of terms in the partial sums. For the Drell-Yan exponent, the addition of perturbative terms 
up to the optimum number N(t( ft)) makes the exponent progressively scale-invariant. 

Smooth behavior under scale transformations is not a property of the resummation of 
Ref . [|[ . If an IR cutoff is used, conventional RG-scale invariance looses much of its meaning 
since an additional and arbitrary scale is introduced. We return to this issue in section 5 
when we discuss the physical cross section and compare scale variations in our approach with 
the corresponding ones of Ref. . 
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4 Leading and subleading logarithmic structures 



In this section we describe the resummed partonic cross section in momentum space, derive 
a perturbative approximation in a specific region of momentum space, and discuss in general 
how to include leading logarithmic structures consistently in the resummation. We also 
demonstrate how the constant gW is determined from the one-loop calculation, and we 
discuss the issue of universality of the leading logarithms. In particular, we demonstrate 
that retention of subleading pieces in the exponent, of a structure similar to the Drell-Yan 
case, would not account correctly for the subleading structure at the one-loop level in top 
quark production. We justify the truncations we use consistently with universality of initial- 
state radiation. 

The general form of the PVR cross section in the MS scheme is 



1 

l-4(l+r?)+VI+^ 



a^ v {r],m 2 )= I dzH(z,a)a' i:j (r],z,m 2 



where the kernel of the hard part is 

H(z, a) = 1 + / dxe E{x ' a) £ Q,{x, a), (89) 

J ■ C\ 



3=0 



and 

x, = In I 

\-z, 

The functions Qj(x,a) in Eq. (|89D appear during the inversion of the Mellin transform 



= In (— !— ) . (90) 



n «-> z 



IT| . They are obtained from the generating function 

f 1 00 pmi pm 2 r> m N+i 

S[P lj P 2l ...P, +1 ]ERei E ^_^±Y 

1 » lima ,3^ JV+1 =o m 1 \m 2 \...m N+1 \ 



/ o \ mi+2m2+...(N+l)m N+ i 

xlim^o(^) -'"-I'd +<-)} (9.1) 



through the identification 



Q[P ll P 2 ,...P N+1 ]=J2Qj ■ (92) 

j=0 
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This notation is slightly different from that of Ref . []TT] . In particular, in the present notation 
the real part is extracted explicitly from the generating function. In Eqs. (|9T| ) and (^) , 



P k = P k {x,a) 



d k E(x,a) 
k\dx k 



(93) 



The function Qj = Qj(x,a) is defined as the set of terms in Eq. ([89]) that contribute j 
more powers of a than of x. The structures Qj resum the j-th subleading logarithms in the 
physical cross section through the index identity 



that connects Eqs. ([Jl]) and fl92|). 

4.1 The leading resummed perturbative cross section and its range 
of validity 

For best accuracy, in a process like the Drell-Yan process, the expression for the hard part 
should include enough of the set of Q/s to reproduce the resummable finite-order subleading 
logarithms up to two-loops. On the other hand, owing to the constrained universality that 
characterizes the ti cross section, we should include only terms that do not create non- 
universal subleading logarithms. In addition, we must discuss the range in momentum space 
where our perturbative resummation is valid. 

We turn first to the issue of the perturbative regime in momentum space. Specification 
of this region follows from general expressions for the inversion of the Mellin transform and 
the meaning of the successive terms in this inversion, once their perturbative approximations 
are used. The kernel of the hard part is provided in Eq. (|89"D , and the functions Qq and Q±, 
obtained from Eqs. (|9~l|) and (|9~3"D, are 




(94) 



k=2 




(95) 
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and 




(96) 
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The functions * = ^ (0) and are the usual Polygamma functions. For simplicity we 
include in the expression for Qi only terms that generate corrections starting at 0(a). 
According to the general discussion at the beginning this section, Q\ contributes one less 
power of x than of a in the integrand of Eq. (|8^), and it is formally subleading relative to 
the contribution of Q . Nevertheless, from Eqs. fl9"o] ) and fl9"6|), we see that this suppression 
is not true for values of x such that P\(x, a) ~ 1. When sin(7rP 1 ) ~ 0, the dominance of Qq 
over Qi is destroyed. For values of x such that P\ > 1, Qq and Q\ are out of phase, with 
amplitudes that are not constrained, and the perturbative dominance of Q over Q 1 is again 
vitiated. 

We conclude that the perturbative region in momentum space is defined by the inequality 
constraint 

P 1 (x z ,a)<l. (97) 

In our discussion in moment space in section 3, we saw that the perturbative approximation 
for the exponent is valid in the interval 1 < Inn < t — l/2a&2, i-e., where terms containing 
equal powers of Inn and a are at most 0(1). Our Eq. fl9~7|) translates this condition con- 
sistently into momentum space. In Fig. 5 the upper two curves show the range in Inn = x 
where the constraint of Eq. ( |9"TD is satisfied. The lower three curves in Fig. 5 pertain to the 
corresponding constraint for the hard part in moment space, Eq. ([31]). We see that the two 
criteria are in good agreement: the value of Inn = x at which Eq. ( |9~T[ ) is satisfied agrees 
fairly well with the value of n at which Eq. (|8l|) is also. 

We take up next the issue of universal logarithmic structures in our resummation. Power- 
counting, discussed generically in section 2.2, must be altered when applied to the kernel 
of the integral in Eq. (^) because the extra integration supplies one more power of the 
logarithm. We begin with the following expression for the kernel of the hard part: 

H(z, a) ~ 1 + r dxe E{x > a) Q (x } a) , (98) 
Jo 

and we rewrite 

Qo(x, a) = Qo(x, a) — Pi(x, a) + Pi(x, a) = A(x, a) + Pi(x, a) . (99) 
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Using Eq. fl63|), we can express 



oo 2m 

e E(x,a) = Y^ a m £(l,m)x l , (100) 

m=0 l=m+l 

and, from the definitions of Qq and Pi, we find the expansion 

oo 

A(x,a) = 5>"M m x m (101) 



m=2 



Therefore, 



/ °dxe E ^A(x,a) = £ a m+2 J2 £ 

m=0 k=0l=k+ 



m 2fe e(Z, k)S, 



fc+1 3 + m + / - k 



xi +m+l ~ k . (102) 



For a fixed m, the maximum monomial in powers of x z is a 2+m x 2m+3 . In the notation of 
section 2.2, the integral in Eq. (|102[) contributes at most to the first-subleading logarithms 



c(2m — l,m). Since we intend to resum leading logarithms only, we ignore the contribution 
from the integral Eq. (|102p . For the same reason, we discard all Q/s but Qo in Eq. (|89f) . 

Our main result for the perturbative resummed partonic cross section, denoted by afj rt , 
is therefore 



af" rt {r],m 2 ) = / dzl+ dxe E( - x ^ P 1 (x, a) 

Jl-4(l+» ? )+4 v / T+^ L JO 



0-^.(7/, z,m 2 ) 



/ dze E(x *> a) aUri,z,m 2 ) . (103) 

il-4(l+7?)+4^I+^ 



The upper limit of integration z = z Q (a) is determined from the equation for the perturba- 
tive regime 

' e* , P 1 (x ,a) = l. (104) 



l-zo 

The preceding expressions in this subsection are valid in either factorization scheme. The 
only distinction is the value of the multiplicative constant in Eq. (|64]). This constant 
can be determined from a comparison of the 0(a) expansion of the resummed partonic cross 
section with the finite-order calculation of Refs. [1], . 
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It is useful for purposes of comparison to expand the cross section to 0(a 2 ). In the MS 
factorization scheme, using Eqs. and (|64|), we obtain 

(105) 



£ [2] (x, a) = ag^x 2 + 2a 2 g { % 2 x 3 /3 . 



Substituting Eq. (|105|) into Eq. ( 103 ), we find that the corresponding finite-order part of the 
resummed cross section is 



o\\r],m 2 ) 



dzll + ag {1) hi 2 f-^- 

pert "'l-4(l+r ? )+4 v / T+^ I V 1 — Z 



1-2 



3 V 1 — z 



(106) 



To determine g^ l \ one can use the corresponding O(a) expression for the cross section 
from Ref. ||, or derive an asymptotic formula near threshold (77 — > 0) for the 0(a) piece of 
the integration Eq. (|106| ) and compare with the corresponding explicit expression of Ref. |L|. 
The former is simpler for the present purposes. In section 4.2, we derive the asymptotic 
threshold formula to demonstrate universality issues. The result is 



.'/.'/ 



(107) 



the same as in the Drell-Yan case. Our Eq. ( |106|) , including the 0(a 2 ) term, is identical to 
that of Ref. §. 

Similarly, in the DIS factorization scheme, using Eqs. (^) and (|64"1) , we derive 



E [2] (x, a) = ag (1) x 2 /2 + a 2 g {1) b 2 x 3 /2 . 
Hence the cross section to 0(a 2 ) becomes 



(108) 



a\Mr],m 2 ) 



dzl 1 + a- — In 2 ^ 



pert Jl-4(l+»7)+4Vl+r? 

1 



2\ 2 J \l-z 



2 

2 y \l-z 



1 - z 



>a' ij (r ] ,z,m 2 ) 



(109) 



By comparison with the corresponding 0(a) piece of Ref. 0, we obtain 

Ji) - or - n W 



(110) 
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again in agreement with the Drell-Yan case. As before, Eq. ( |109| ), including the 0(a 2 )-term, 
is identical to that of Ref . [Q . 

We close this subsection with simplified analytical expressions for the perturbative regime 
of the cross section, Eq. (|97D . For this purpose, we may use the convergent expressions, 



Eqs. (gg) and (|68Q . In the MS factorization scheme, we obtain 

,(i) 



P x (z a , a) ~ - 9 — ln(l - 2ab 2 x z ) < 1, (111) 



and in the DIS scheme, 



o(i) o(i) 
Fife, a) ~ - V- ln(l - 2a& 2 x z ) + ±— ln(l - afe 2 x 2 ) < 1. (112) 
o 2 o 2 

To obtain transparent analytical expressions, we specialize here to the one-loop approxima- 
tion 

2ab 2 = ln _1 (m/A). (113) 

We obtain 

. a x l-exp(-6 2 /s«) 

l-z>(— ) MS, (114) 

x _ oxp(-(, 2 / 9 ( 1 ))/2 

1 - p( - 62/sllJ)/2 DIS. (115) 

We observe that the non-perturbative regime is suppressed mainly by A/m to the first power, 
with a correction that depends on the color factors of the partonic production process. For 
m = 175 GeV, 

A/m~10~ 3 . (116) 

Because 6 2 is positive (6 2 ~ 2), the non-perturbative regime is an increasing function of 
the color factors. If the partonic cross sections had no color-factor enhancements, < 
1, Eqs. ( 114 ) and (|115| ) would be in close agreement with the perturbative regime of the 



exponent in moment space n < m/A, Eq. (0), with the direct substitution 1/n <-> 1 — z. 
For example, for g^> = 1, 

1 - z > 2 x 10~ 3 MS and 1 - z > 1.5 x 1(T 3 DIS , (117) 
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in good agreement with Eq. ( |116| ). In reality, = 8/3 in the qq channel, therefore 

1-2>2x10" 2 MS and 1 - z > 8 x 1(T 3 DIS . (118) 

For the gg channel, g^ = 6, and we find 

1 - z > 1 x KT 1 MS and 1 - z > 5 x KT 2 DIS . (119) 

These regions are narrower than the estimate based on the perturbative properties of the 
exponent only, especially for the gg channel. The soft gluon region is probed more deeply 
in the DIS scheme than in the MS scheme, closer to A/m. We discuss these properties 
numerically in section 5. 

4.2 Non- universal sub leading logarithms 

In this subsection, we address universality of the logarithmic structures in the threshold 
region more explicitly. It is sufficient to examine the issue in the MS scheme. For complete- 
ness, and to compare with the analytical results of Ref. 0, we derive the near-threshold 
asymptotic properties of the 0(a) term of Eq. ( |106| ). We demonstrate that if we were to 
keep the subleading logarithm in the Drell-Yan exponent, Eqs. fl53|) and (p5|), the resulting 
term in the 0(a) expansion of Eq. (j39| ) would not be the same as the one obtained in the ex- 
plicit next-to-leading order calculation of Ref. 0. This is a demonstration of the consistency 
requirements of the series of truncations that resulted in Eqs. fl63|) , (0) , and (|103|) . 
In the MS scheme, Eq. ([53]), the full 0{a) Drell-Yan exponent is 

eW{x, a) = ag {1 \s 2}1 x 2 + s 1A x + s ,i) , (120) 

where the coefficients are found in Eq. (|55|). Suppose we allow an undetermined linear term, 
anx, that would represent the deviation of this exponent from the Drell-Yan form. We 
discard the constant term in this exercise, because it does not contribute logarithmic terms 
to the cross section at first order. Substituting from Eq. fl5"5|) we have 

E [1] (x, a) ~ a{g {l) x 2 + 2g w ^x + kx), (121) 

where 7 is the Euler-Mascheroni constant. The kernel of the hard part, containing subleading 
logarithmic structures, is given generally by Eq. (^). To account for subleading logarithms 
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with full accuracy in the Drell-Yan case, we must keep Qo and Qi, Eqs. fl9"5|) and 
Expanding to 0(a) and substituting into Eq. (|88|), we obtain 



oj/ (77,777 ) = a 

J Jl-4(l+77)+4v^+»7 



<7« In 2 (-L-) + «ln (^)] 4(r^,m 2 ) . (122) 



The piece of the Drell-Yan expression linear in the logarithm has dropped out. In other 
words, the exact evaluation of the cross section of Eq. ( |122| ) with the Drell-Yan resummation 
expression gives the same answer as the truncated expression, Eq. ( |106| ). This is something 
well-known in the Drell-Yan case and is expressed by absence of the function #2 in the MS 



scheme, Eq. ([T6|). 

The integral representation of the partonic cross section does not contain the subleading 
logarithms of the Drell-Yan case. This is not true for its asymptotic evaluation near thresh- 
old, as we shall see below, owing to the form of the Born cross section. We focus on the qq 
channel only. From Eqs. (|2H|) and (|122|), we obtain 



4 } (77,m 2 ) = -vr 3 « 3 - / dz[l - (1 - z)r\ x 

9 S Jl-4(l+ ?? )+4 v / T+J7 



[1 - (1 -z)t] 2 -At 



2t 



1 - z 



^[1-{1-z)tY-At 
1 



X 



1 - z 



(123) 



where r = m 2 / s. Introducing the variables x = [1 — (1 — z)t] 2 — 4r and[Q f3 2 = 1 — 4r, we 
derive 



a S = 7T a l / dx 
9s Jo 



9 {1) 1^ (t^ 1 " V 1 " ^ ~ X)] ) +Kln (l^ [1 " V 1 " ^ " *)]) 



.(124) 



The threshold is exhibited at (3 — > 0. Expanding 



1 - yjl - (/3 2 - x) ~ {(3 2 - x)/2 



(125) 
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and using 



ln fc x = lim f 7T- | x e , 
e-+o I de 



(126) 



we find 



(i) 4 3 /? 3 

qq — n s 9 

9 mr 



d 



5(3/2, 1 + e) + 1—^5(1/2, 1 + e) 



2/3 2 



9 / 2(3 2 



de \ l-(3 2 



5(3/2, 1 + £) + L^B(l/2, 1 + e) 



(127) 



After computing the derivatives, we obtain 



w " 9 s m 2 ^ 



I In 2 ( 7^ ( 5(3/2, 1) + ^5(1/2, 1) 



l-(3 2 ) \ 



2p 2 



+21n (y^) (b(3/2, 1)(¥(1) - ¥(5/2)) + ^^(1/2, 1)(¥(1) - ¥(3/2)/ 
+5(3/2, 1) (($(1) - ¥(5/2)) 2 + ¥'(1) - $'(5/2)) 



2/3 2 



5(1/2, 1) (¥(1) - ¥(3/2)) 2 + ¥'(1) - ¥'(3/2) 



1" I t^L^I ( 5(3/2, 1) + 1—^5(1/2, 1) 



l-(3 2 J \ 



+5(3/2)(¥(l)-¥(5/2)) 



2/5 2 



2/5 2 



5(1/2, !)(¥(!)- ¥(3/2)) 



(128) 



Discarding the non-logarithmic pieces and using the relations ¥(1) = —7, ¥(3/2) = —7 — 
In 4 + 1, we observe that all the coefficients of the logarithms can be rational. Introducing JIJ 
p = At, we end with 



v. 



(1) _ a[^r 
qq ~ m 2 9 1 



ln 2 (8/5 2 ) - [2gW-K] ln(8/3 ; 



(129) 



Our Eq. ( |129| ) may be compared with the result of Ref. |]J in the MS scheme, denoted a^f, 



3/(1). 



99 • 



" m 2 9 13 



41 1 
ln 2 (8/3 2 )--ln(8/3 2 )] 



(130) 
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We conclude that, for the leading logarithm, 

g W = 8/3 = 2C F , (131) 

the same value found in Eq. ( |107| ). However, for the next-to- leading linear logarithm, 

2g (1) = 16/3 ^ 41/6 . (132) 

A non-zero value for k would be required. This exercise shows that universality between the 
Drell-Yan case and ti production is restricted to the leading logarithms only, as must be the 
resummation of these structures. The subleading logarithms in the full exponent of Eq. ([53]) 
are uncertain. 

The analysis of this subsection justifies our truncated expressions, and it motivates sev- 
eral equivalent ways to define a quantitative measure of the overall theoretical uncertainty. 
The most straightforward is to vary the renormalization/factorization (hard) scale within 
a reasonable range and to use the variation of the cross section as a measure of the non- 
universal subleading structures. It is important to recall, however, that truncation of the 
resummation of leading logarithms at a given higher order produces subleading logarithms 
at that order, as the 0(a 2 ) terms in Eqs. ( |106|) and (|109|) show. We will assume that these 



particular subleading logarithms are approximately universal because they come from the 
evolution properties of the resummed universal leading logarithms of lower orders. This is 
true for the MS Drell-Yan cross section at two loops (16| . 



To summarize, the resummation exponent in the perturbative regime in the MS factor- 
ization scheme is 

2CV,| £ a^^^x^-Hm/fi) £ a^f-^—x^ . (133) 

In the DIS factorization scheme, it is 
E(x,a(fi),fi,N(t(iJ,))) = 

( N(t(fM))+i h p ~ lc ?p w(to))+i hf~hp ) 

2C y | E ^)-^Yf +l -Hrn/^) E ^P-j-^ 
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(N{2t(jM))+l h p-l N(2t{fi))+1 h p-l n 

-2CJ £ ^)-^—-x^-\n(m/^ £ ocO{h)°^x» . (134) 
I pi P(P + 1) p P J 

The upper integer in the two sums is N(t(n)) = [t(/x) — 3/2], and t(fi) = l/(2a(/i)6 2 )- 
The perturbative resummed cross section in the corresponding scheme is obtained from the 
general expression Eq. fl8U|), within the region defined by Eq. (0). This restriction can be 
incorporated if we define zq through 

P 1 fin ( - i —) , a(ji),m/n, N(t(fi))) = 1 . (135) 

V \1- Zo{fji,m/fji)J J 

Imposing the constraint that only leading logarithms are resummed, we arrive at our per- 
turbative resummed cross section 

< rt (^™ 2 ,/! 2 ) = r (P ' mM _dze E ^ a ^ m ^ N ^aU V ,z,m^ . (136) 

Jl-4(l+»y)+4 v / l+r ? 

The derivatives of the Born cross sections, a^r], z,m 2 ), are found in Eqs. (|20| ) and (p2j). 

We use the formulas above to calculate our predictions for the perturbative resummed 
cross section. The rest of the phase space in Eq. ( |136| ), from z = zq to the threshold 
z — 1, is non-perturbative, and we declare ignorance about the cross section in that region. 
(In section 6 we engage in some speculations on non-perturbative physics.) A measure of 
the dominance of the universal logarithmic structures and their resummation in the per- 
turbative regime in our approach is provided by the change in the cross section when the 
renormalization/factorization scale is varied within a reasonable range around the "central" 
value fi = m. This variation measures the importance of subleading logarithmic structures 
that are non-universal and non-resummable in this approach. We use variation with n as a 
working hypothesis for the bulk of the perturbative uncertainty^]. 

5 The resummed cross section 

In this section we present analytical and numerical results for the resummed cross sections. 
Section 5.1 is devoted to the partonic cross sections as functions of the variable r], while in 
section 5.2 we present the physical cross sections obtained after convolution of the partonic 
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cross sections with parton densities. We furnish details on the technicalities involved, both 
kinematical and numerical. In section 5.3 we make comparisons with the predictions of 
Laenen, Smith and van Neerven|3|, [3J], denoted LSvN. 



5.1 Partonic cross sections 

The resummed partonic cross section has the form 



< rt (r/,m 2 ,/i 2 ) = r° (M,m/M) ^^W^Woj^^^), (137) 

Jl-4(l+ ?? )+Vl+r ? J 



The upper limit of integration, z (/i, m//i), is provided by Eq. ( |135| ). For physical but rel- 
atively large values of 77, the lower limit of integration may become negative, unlike the 
Drell-Yan process in which < z < 1. In this situation, far away from threshold, resumma- 
tion of initial-state gluon radiation is irrelevant, and we do not perform resummation outside 
the range < z < 1. Finding the roots of the lower limit of integration in Eq. ( |137| ) we see 
that this equation is unconstrained for 



rj < (1 + V2)74 - 1 = 0.4571068 . 



138) 



Equation (|138 ) defines the region in phase space in which threshold effects are important. 
Above this value of 77, Eq. (|137|) should be constrained so that it includes only the phase 
space where gluon radiation is produced near threshold. 

In order to achieve the best accuracy available we wish to include in our predictions as 
much as is known theoretically. The exact next-to-leading order (one-loop) cross section 
is known Jl], |2j, but the full two-loop calculation does not exist. Thus the best we can do 
at present is to include the full content of the one-loop partonic cross section along with 
our resummation of the leading logarithms to all orders. In so doing, we include both 
non-universal one-loop subleading logarithms and constants. This procedure is in common 
with previous resummations of this process^] and the Drell-Yan pro cess |12]. Our "final" 
resummed partonic cross section can therefore be written [H 



pert/ 2 2 

afj (77, m ,// 



a. 



pert 



a (0+l)^ ...2 ,,2, 



(77, m ,// 



+ 4 0+1) (r?,m 2 , /U 2 ) . (139) 



pert 
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The second term is the part of the partonic cross section up to one-loop that is included in 
the resummation, while the last term is the exact one-loop cross section]!], |2|. 

For numerical purposes it is best to convert the "improper" integrations induced in 
Eq. ( |137p by the form of the Born cross sections, cr^(r/, z,m 2 ) (see Eqs. (|20|) and fl22|)), to 

integrations without a numerical (but analytically integrable) singularity. This is achieved 
by the transformation 



* = f- p- T ft -,)). - < 140 > 

In addition, it is easier to parametrize the perturbative regime, Zo, by the moment variable 
n : 

—!— = no, n = e X0 , (141) 
l — z 

where x is the first root of the equation 

P 1 (x (fi,m/tJ>),a(fM) 1 m/fi 1 N(t(ji))) = 1 . (142) 

Taking into account the constraint on rj, Eq. ( |138| ), we can write the various pieces of 
Eq. ( P9D as follows: 

fU (r.no) 

af/ rt (^,m 2 ,/i 2 ) = / dxe E ^' a ^ m ^' N ^^a i:j {x) , (143) 

JL(t) 

where 



U(r, noC/i, m/n)) = ^1 - (1 _ r/ ^ m//l))2 , (144) 



L(r)=e( T I x Wl- - 4r NO , (145) 

W Ul + v/2) 2 / V (l-^) 2 



and 



r\/l — x 2 



gW = ln i ^-2^ J • (146> 

The second term in Eq. ( |139| ) is obtained from Eq. (|143| ) by expanding the exponential up 
to 0(a). The exponents for the two factorization schemes are found in Eqs. ( |133j ) and ( |134| ). 
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while the transformed functions &ij(x) are 



. . 32 3 2 r 3 / 2 / 1+x 2 



(1 -x 2 ) 5 ^ 



(147) 



and 



1 r 1 / 2 

~ / \ 3 2 

°gg\ x ) = 2 n a 



S (l-X 2 ) 3 /2 



-x 2 + 5r 



. ( 2x 4 - -a; 2 + 3rx 2 + 2 - 3r)xln(^^- 
3 / \ 1 — x 



| --x 4 - 4x 2 + tx 2 + - - 3r 



12r 

1 -x 2 



148) 



In Fig. 6 we show the perturbative boundary no(fi,m/fjL), obtained from Eqs. ( 141|) and 
( |142| ) , as a function of the hard scale /i, for m = 175 GeV. The scale /z is an artifact of pertur- 
bation theory, and /^-variation is associated with truncation of the perturbative expansion. 
In a resummation such as ours, all significant perturbative knowledge of threshold effects is 
exhausted (all large perturbative threshold corrections are included), and the perturbative 
regime we calculate should also be insensitive to artifacts such as /i. Therefore, one expects 
the perturbative boundary to be independent of the hard scale /i and to depend only on 
the physical scale characterizing the threshold, m. We observe that the full Drell-Yan ex- 
ponent produces a function Pi that is almost exactly hard-scale invariant, establishing that 
our resummation conforms to this intuition. On the other hand, the truncation to universal 
logarithmic structures that we use for our predictions shows some scale dependence. 

In Fig. 7 we show the resummed partonic cross sections afj Tt as a function of rj, for 

\i = m and m = 175 GeV, for both production channels in the MS factorization scheme. 
We choose a logarithmic scale in rj to expand the threshold region. We also show the lowest 
order and next-to-leading order counterparts. The three curves differ substantially in the 
partonic threshold region r\ < 1, with the final resummed curve exceeding the other two. 
Above rj ~ 1, our resummed cross sections are essentially identical to the next-to-leading 
order cross sections, as is to be expected since the near-threshold enhancements that concern 
us in this paper are not relevant at large rj. In both the qq and the gg channels, we note that 
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the size of the 0(a^) term exceeds that of the 0(a 2 s ) term for 77 ~ 0.1, and the ratio grows 
as 7] decreases. This behavior is contrary to the notion underlying perturbation theory, that 
successive terms in the perturbation series should be smaller, and is cited in our Introduction 
as the motivation for resummation at small 77. 

It is useful to translate our definition of the perturbative regime directly into a statement 
about the perturbative region in 77. Our perturbative resummation probes the threshold 
down to the point 

V > Vo = (149) 
2n 



where no is calculated from Eqs. ( |141| ) and ( 142 ). Below this value, perturbation theory 



resummed or otherwise, is not to be trusted. For our central predictions we choose to accept 
the exact next-to-leading order results throughout the phase space, but the non-perturbative 
region is a source of non-perturbative uncertainty, subject to model-building. At m — 175 
GeV, our resummed cross sections become identical to the next-to-leading order cross sections 
below T) > 0.008 for the qq channel and 77 > 0.05 for the gg channel, a consequence of our 
decision to restrict resummation to the perturbative domain. The difference reflects the 
larger color factor in the gg case. 

We return to possible non-perturbative contributions in section 6. 

5.2 The physical cross sections 

We use S to denote the square of the hadronic center-of-mass energy. Once the partonic cross 
sections are calculated, the physical cross section for each production channel is obtained 
through the factorization theorem. 



Am 2 r-r^r- 1 
cr i:j (S,m ) = —g- dri<& 



'j 



Am 



2 



S 

The parton flux is a convolution of parton distributions 



1 + 7?U 5 



Oij (77,771 ,fi )• (150) 



®ij[y^ 2 }= —fi/hA x ^ 2 )fj/h 2 (y/x^ 2 ). (151) 

Jy X 



We use CTEQ3 parton distributions |15| in the appropriate factorization scheme. The total 
physical cross section is obtained after incoherent addition of the contributions from the the 
qq and gg production channels. We ignore the small contribution from the qg channel. 
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A quantity of phenomenological interest is the differential cross section 



daij(S,m 2 ,ri) 4m 2 

Tr) = ~s~ ij 



Am 
~S 



2 



[i + v), V 



0-^(77, m 2 ,/i 2 ) . (152) 



The differential distribution is a RG-invariant quantity, and it is perhaps measurable. Its 
integral over rj is, of course, the total cross section. In Fig. 8 we plot these distributions for 

the two production channels for m = 175 GeV, v^S* = 1.8 TeV and fi = m. Convolution with 
the parton flux enhances the relative importance of the region of small rj. We observe that, 
at the energy of the Tevatron, resummation is significant for the qq channel and less so for 
the gg channel. Figure 8 is useful for back-of-the-envelope estimates of possible contributions 
from the non-perturbative regime, discussed in section 6. 

We show the total tt-production cross section for various energies in Fig. 9, and in Table 1 
we provide numerical values. The central value of our predictions is obtained with the choice 
fj,/m = 1, and the lower and upper limits are our estimate of the perturbative uncertainty. 
These upper and lower values are the maximum and minimum of the cross section in the 
range of the hard scale \xjm G {0.5, 2}. For the range of top quark mass shown, the minimum 
occurs at fi/m = 2, while the maximum occurs at fi/m ~ 0.7, as is also shown in Fig. 9b. 



Our prediction of Fig. 9a is in reasonable agreement with the published data|T7|]. The 
/^-variation of our resummed cross section, shown in Fig. 9b is smaller than that of the 
next-to-leading order cross section, as expected for an all-orders resummation. (Owing to a 
computer compiler error, the variation with /i of the next-to-leading order results in the first 
paper of Ref. || is incorrect. Our overall resummed predictions are essentially unchanged. 
The variation shown in Fig. 9b should be correct.) 

Our resummed cross sections at VS = 1.8 TeV are about 9% above their next-to-leading 
order counterparts computed with the same parton distributions. To gain numerical insight 
into the magnitude of this increase, we may examine the growth of the cross section in the 
dominant qq channel that would be expected in a series of fixed-order calculations if only 
leading- logartihmic threshold contributions are included in successive orders in a. We take 
the universal leading- logarithmic contributions at 0(a) above the Born level from the full 
next-to-leading order calculations of ti production and the leading-logarithmic contributions 
at 0(a 2 ) from calculations of the Drell-Yan cross section. For this exercise only, we set aside 
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our perturbative constraint, Eq. ( |135| ), and instead we integrate over the entire threshold 
region < z < 1. We present the resulting finite order physical cross sections in the MS 
scheme in Fig. 10. (For the Born cross section, we integrate over all z.) We observe that 
the cross section 0-(° +1 ) obtained from the leading- logarithm terms, only, through 0(a) is 
in remarkable agreement with the cross section obtained from our full next-to-leading order 
calculation for ti production. Moreover, our predicted resummed cross section lies part way 
between o"^ 0+1 - ) and the cross section o-( 0+1 + 2 ) obtained from the leading-logarithm terms, 
only, through O (a 2 ). At m = 175 GeV, the increase of a^ 04 ^ over the Born result is a 22% 
effect, and the further increase of cr(°+ 1+2 ) over <j( 0+1 ) is another 14% effect. We conclude 
that the roughly 9% increase of our final resummed cross section above the next-to-leading 
order cross section is quite reasonable. It would be surprising if it were much less. 

Table 1: The total ti production cross section at y/S = 1.8 TeV and its perturbative un- 
certainty. The theoretical error band, normalized with respect to the central value, represents 
an almost constant uncertainty of 9 to 10%. 



m (GeV) 


tfert ( min ; P b ) 


a pert (central; pb) 


<eri ( maX ; P b ) 


150 


11.76 


12.72 


12.90 


155 


9.87 


10.68 


10.83 


160 


8.33 


9.01 


9.14 


165 


7.06 


7.63 


7.73 


170 


6.00 


6.48 


6.57 


175 


5.10 


5.52 


5.59 


180 


4.36 


4.71 


4.78 


185 


3.73 


4.04 


4.09 


190 


3.20 


3.46 


3.51 


195 


2.75 


2.98 


3.02 


200 


2.37 


2.57 


2.60 


205 


2.04 


2.21 


2.24 


210 


1.77 


1.91 


1.94 


215 


1.53 


1.65 


1.68 


220 


1.32 


1.43 


1.45 


225 


1.15 


1.24 


1.26 


230 


0.99 


1.08 


1.10 


235 


0.86 


0.94 


0.96 


240 


0.75 


0.81 


0.83 


245 


0.65 


0.71 


0.72 


250 


0.57 


0.62 


0.63 
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As remarked earlier in the paper, our resummation includes only the universal logarithmic 
structures. It is reasonable to inquire whether and by how much the predicted cross section 
would change if subleading logarithmic structures are included. One good representation of 



the possible effect of subleading structures is the full Drell-Yan exponent itself, Eq. (|86|) . As 
shown in Fig. 4, this exponent is larger than our leading exponent, but the corresponding 
perturbative regime, calculated through Eq. ( |135| ), is smaller. To obtain a cross section that 
includes the Drell-Yan subleading logarithms we use Eqs. (p8|)-(|90|) with the sum over Q/s 
replaced by the full first term, Qq. It is shown as the dashed line in Fig. 11 for the qq channel. 
We see that the curves differ little from each other, and that the Drell-Yan resummation 
prediction is within our uncertainty band. This example substantiates our belief that /i- 
variation is an adequate measure of perturbative uncertainty and that it includes the effect 
of non-resummable, non-universal logarithms. 

A second source of uncertainty, that is partly phenomenological, partly perturbative, and 
partly correlated with non-universal logarithms, is associated with the use of different parton 
distributions. The parton set we use is a next-to-leading order determination of the quark 
and gluon densities. Since we use resummed partonic cross sections, it is arguably true that 
we should also use parton densities based on resummed expressions for deep inelastic lep- 
ton scattering cross sections and other processes used in the determination of the densities. 
However, no such densities exist. It is common practice phenomenologically to repeat calcu- 
lations with different sets of parton and to estimate thereby a second source of uncertainty. 
Except for the fact that different data sets are used, or are emphasized differently, in different 
determinations of parton densities, and that the fitting programs differ, we opine that the 
practice of adding in quadrature, or otherwise, uncertainties associated with fi variation and 
those associated with different parton sets involves significant double counting. Much of the 
difference among modern parton sets reduces to a difference in A which largely affects a(fi). 
Thus, this difference is correlated with the /i-variation we consider above, and should not be 
treated as an independent error. Over the range fi/m G {0.5,2}, the band of variation of 
the strong coupling strength a s is a generous ±10% at m — 175 GeV. 

In Fig. 9c we present our predictions for an upgraded Tevatron operating at \/S = 2 
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TeV. Our cross section is larger than the next-to- leading order one by about 9%. We predict 

a tl {m = 175 GeV, VS = 2 TeV) = 7.56^ pb . (153) 

At m = 175 GeV, the value of the cross section at = 2 TeV is about 37% greater than 

that at y/S = 1.8 TeV. 

Turning to pp scattering at the energies of the Large Hadron Collider (LHC) at CERN, we 
note a few significant differences from pp scattering at the energy of the Fermilab Tevatron. 
The dominance of the qq production channel at the Tevatron is replaced by gg dominance 

at the LHC. Owing to the much larger value of a/S*, the near-threshold region in the suben- 
ergy variable is relatively less important, reducing the significance of initial-state soft gluon 
radiation. Lastly, physics in the region of large partonic subenergy ^/s, where straightfor- 
ward next-to-leading order QCD is also inadequate [|I8||, may become more significant for 
ti production at LHC energies than the effects of initial-state radiation. The approach of 
our paper is limited to the resummation of initial-state gluon radiation only. We present 
estimates in Fig. 9d of the cross section for LHC energies of 10 and 14 TeV. We obtain 

a a (m = 175 GeV, Vs = 14 TeV) = 760 pb . (154) 
5.3 Comparisons with other calculations 

In earlier sections we comment on differences between our formalism and that of Laenen, 
Smith, and vanNeerven (LSvN) |3|, [|. Here we compare aspects of our numerical predictions. 
The comparison most relevant to experiment is that of our Table 1 and the corresponding 
table in Ref. ||. Our central values are 10 to 14% larger (the difference increases with 
mass), and our estimated theoretical uncertainty is 9 to 10% compared with their 28% to 
20% (decreasing with mass). The two predictions have overlapping uncertainties and are, 
in this sense, in agreement. In commenting on differences, we remark that our Born cross 
section is about 3 to 5% larger than LSvN's Born cross section. The difference arises from 
the different parton distributions used in the two calculations, including differences in A 
which alone account for half or more of this increase. This source of difference should be 
kept in mind when comparisons are made with next-to-leading order calculations and various 
previous results. 
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It is important to stress that the theoretical uncertainties are estimated in quite different 
ways in the two methods. We use the standard /i-variation, whereas LSvN obtain their 
uncertainty primarily from variations of their undetermined IR cutoffs. From a theoretical 
point of view, a study of the variation of the cross section with the hard-scale /i is important 
because it deals with stability of the calculation under variation of a perturbative, but not 
directly determinable renormalization/factorization scale. This statement applies as well to 
the LSvN calculation, above and beyond the choice of their IR cutoff. The role of the IR 
cutoff is to measure ignorance of non-perturbative effects in the LSvN approach. One of 
the advantages of a resummation calculation should be diminished dependence of the cross 
section on /x, less variation than is present in fixed-order calculations. 

It is also instructive to compare the numerical values of our perturbative boundary, 
r]o(Tn, fi/m), with the corresponding cut in rj produced by LSvN's IR cutoff fio- In the MS 
scheme, their cut is0 

v(l*o,v) = \ {~) ' ( 155 ) 

and in the DIS scheme, it is 

V(IM>,IJ) = \ {~) ■ ( 156 ) 

Values of t}(liq, p,) are provided in Tables 2 and 3 and are compared with our perturbative 
boundary r) (m, fx/m). 

In the paragraphs to follow, we compare /i-variation, scheme dependence, and the influ- 
ence of the difference between our perturbative boundary and LSvN's cut in rj. In Tables 2 
and 3, we reproduce numbers from Ref. H and show the corresponding values we calculate. 
The comparison is made at m =100 GeV in Table 2 because it is at this value of the top 
mass that they provide results in the MS scheme for both the qq and gg channels. 
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Table 2: Physical cross sections in pb for m = 100 GeV at VS = 1.8 TeV. The LSvN 
predictions are shown for various choices of their IR cutoff /x - Our perturbative boundary 
rjo and LSvN's phase space cut 77(^0, (J>) are also shown. Absences in the LSvN entries denote 
extremely large numbers. 



<j qq (m = 100 GeV; MS) 


fi/m = 0.5 


fj,/m = 1 


\xjm = 2 


Born 


64.3 


48.4 


37.3 


NLO 


54.6 


55.0 


51.3 


pert 


56.9 


59.4 


57.2 


LSvN(/i = 0.1m) 


52.6 


88.9 




LSvN(/i = 0.2m) 


42.5 


66.4 


110.5 




0.007 


0.011 


0.014 


r](/i = 0.1m, 11) 


0.004 


0.0005 


0.000063 


77(^0 = 0.2m, fj) 


0.032 


0.004 


0.0005 


a ( m = 100 GeV; MS) 


\xjra = 0.5 


\xjm = 1 


\ijm = 2 


Born 


36.4 


23.8 


16.1 


NLO 


44.5 


40.5 


34.0 


pert 


44.8 


41.6 


35.0 


LSvN(/i = 0.2m) 


34.3 


88.7 




LSvN^q = 0.3m) 


28.1 


42.4 


820.0 




0.036 


0.063 


0.1 


77(^0 = 0.2m, n) 


0.032 


0.004 


0.0005 


r)(n = 0.3m, fj) 


0.108 


0.0135 


0.0017 


a tl {m = 100 GeV; MS) 


fi/m = 0.5 


n/m = 1 


n/m = 2 


Born 


100.7 


72.2 


53.4 


NLO 


99.1 


95.5 


87.4 


pert 
a tt 


101.7 


101.0 


92.2 


LSvN^q = {0.2,0.3}m) 


70.6 


108.8 


930.5 
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Table 3: Physical cross sections in pb and perturbative boundaries at m = 150 GeV and 
yS = 1.8 TeV. The corresponding LSvN predictions are also shown. 



qqyiil — LOU Vjrt: V , IJ lO J 


fji 1 ill — yj.o 


ill >y) 1 

11 III 1 


III YYl / 

11 I III Zi 


ATT /~\ 

1NLU 


9.42 


9.31 


8.57 


vert 


9.76 


9.92 


9.31 


T C AT/ C\ 1 \ 

LbvlN(/io = 0.1m) 


7.9 


10.0 


9.7 


%0) 


0.0017 


0.0024 


0.0028 


?7(/io = 0.1m, fx) 


0.02 


0.005 


0.00125 


a gg (m = 150 GeV; MS) 


fx/m = 0.5 


Lx/m = 1 


Lx/m = 2 


NLO 


2.51 


2.22 


1.81 


pert 
99 


2.53 


2.30 


1.89 


LSvN(/i = 0.2m) 


1.76 


4.38 






0.033 


0.055 


0.083 


r)(fx = 0.2m, fx) 


0.032 


0.004 


0.0005 



Tables 2 and 3 show that our resummed cross sections satisfy the test of stability under 
variation of the hard scale fx. The resummed results show less variation than the next- 
to-leading order cross section. On the other hand, this is not true of the resummation of 
Ref. 0. This distinction is linked to the absence of undetermined IR cutoffs in our method 
and the specific RG-invariant exponent we use. Both of these differences contribute to the 
instability apparent in the results of Ref. 0j. The LSvN results show less variation with ix 
in the DIS scheme than in the MS scheme, the reason they provide their final predictions in 
the DIS scheme. As shown in Table 3, the variation with \x of their qq cross section in the 
DIS scheme is about 21%. For comparison, the next-to-leading order cross section shows a 
variation of 9% and our resummed cross section a variation of 6%. 

Our perturbative boundary, i] (n), is fairly insensitive to fx- variation. In the MS scheme, 
it changes by about 30% around its central value (for fx = m) while fx itself changes by 50 to 
100% from the central value. In the DIS scheme the changes in r)o(fx) are even smaller, ranging 
from 16 to 28%. The mild changes in r/o(/x) accord with our earlier physical expectation that 
the boundary is characteristic of the particle's mass, not of the perturbative artifact fx. If a 
complete resummation were possible, as in the Drell-Yan case shown in Fig. 6, our boundary 
r] would be a function of mass only. We may contrast the modest changes of our tfo(fx) 
with the fact that the IR boundary of LSvN varies by three orders of magnitude in the MS 
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scheme. This large change is partly responsible for the unstable behavior of LSvN's cross 
section with //. In the DIS scheme, the IR boundary varies by an order of magnitude for 
fiQ = 0.1m (the value used for their central value prediction in Ref. [f|]), and at /i/m = 1 it 
is twice as large as our r] Q . The fact that these quantities are of the same order of magnitude 
makes LSvN's and our final predictions for the qq cross section comparable. Since the qq 
channel is dominant, our final predictions for the total ti cross section at the Tevatron are 
also equal within uncertainties. 

Scheme dependence is an extra source of theoretical uncertainty, but it should produce 
minimal differences for physical cross sections. The results in Tables 2 and 3 show that 
this is not true of the the LSvN cross sections. We provide our main predictions in the MS 
factorization scheme^]. To check for possible scheme dependent uncertainty, we perform 
our resummation for the dominant qq channel in both schemes. The cross sections presented 
in Table 4 show that scheme dependence is insignificant in our approach, resulting in a 
difference of about 4% for the cross section. 

Table 4. Physical cross sections in pb for the qq channel: DIS versus MS scheme. 



m (GeV) 


<t, ? (DIS; MS ) 


\xjm = 0.5 


yu/m = 1 


fi/m = 2 


100 


NLO 


52.3 


54.6 


52.4 


55.0 


48.9 


51.3 




pert 


53.9 


56.8 


55.5 


59.4 


52.8 


57.2 


125 


NLO 


21.1 


21.9 


21.0 


21.8 


19.5 


20.1 




pert 


21.9 


22.9 


22.3 


23.7 


21.1 


22.6 


150 


NLO 


9.42 


9.68 


9.31 


9.53 


8.57 


8.73 




pert 


9.76; 10.16 


9.92; 10.42 


9.31 


9.87 


175 


NLO 


4.46; 4.54 


4.39; 4.43 


4.01 


4.02 




pert 


4.63; 4.78 


4.69; 4.87 


4.37 


4.58 


200 


NLO 


2.20; 2.21 


2.15; 2.14 


1.96 


1.93 




pert 


2.29; 2.34 


2.30; 2.37 


2.14 


2.21 



In a very recent paper |19| doubts are expressed about the numerical importance of resum- 
mation for top quark production at Fermilab Tevatron energies, and criticisms are leveled at 
our formalism and that of Ref. |J. We consider the criticisms unfounded. The analysis pre- 
sented in our current paper substantiates the work we presented in Ref. || . As demonstrated 
in section 3, our perturbative resummation exponent, e.g., Eq. (|6"5|), contains no factorially 
growing terms in its expansion. The analysis we present of our Figs. 1 through 5 shows that 
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the perturbative region in which we apply resummation remains far removed from the part 
of phase space in which renormalon poles or non-perturbative residual uncertainty could be 
influential. Upon expansion in terms of the QCD coupling strength, fixed at the scale of the 
top quark mass, our formalism produces the universal leading logarithmic structures that 
are found at next-to-leading order (one-loop) in top quark production and at two-loops in 
the Drell-Yan process. This is an all-orders expansion, but it is convergent because its coeffi- 
cients do not have factorial growth, and the momentum scale of these structures is restricted 
to our calculable perturbative regime. Our approach is consistent with, and makes maximum 
use of, the perturbative properties of QCD and established finite-order calculations for the 
relevant processes. With respect to the numerical importance of resummation, we remark 
that the enhancement of the cross section produced by resummation in our calculation is a 

modest 10% above the next-to-leading order result, at a top mass of 175 GeV and v^S* = 1.8 
TeV. That the increase cannot be much less than 10% is suggested simply by an examination 
of the leading-logarithmic contributions at the two-loop level, described in Fig. 10. As shown 
first by LSvNg, it is not necessary to transform to moment space or to resum in order to 
obtain an increase of the cross section in the neighborhood of 10%. 

6 Estimates of Non-Perturbative contributions 

In this section we hazard some estimates of the non-perturbative contributions. This dis- 
cussion is based on educated guesses founded mainly on our intuition of what a physical 
threshold is as well as on the behavior of parton distributions near the threshold. Since 
even the gluon parton distributions are not altogether well established, our estimates may 
be viewed as rough physical understanding expressed quantitatively. 

An examination of Fig. 8 indicates that the non-perturbative regime is very small, es- 
pecially for the qq channel. We focus on this channel first, where the support of the non- 
perturbative interval is 770 ~ 0.008. We recall that the area under the next-to-leading order 
curve is already included in our resummed cross section. Concentrating on the solid curve, 
we assume that the threshold behavior in the non-perturbative regime is a finite continuation 
of the solid curve to r\ = 0. Reasonable guesses include a continuation of the solid curve as 
a constant from its peak value down to rj = 0, or, at most, as a smooth function having the 
same slope as the solid curve just above the peak. In the former scenario we get an extra 
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contribution 

5a q g ~ 5 x 0.008 = 0.04 pb. (157) 
In the latter case we add a further 

5'a q q ~ ^(15 - 10) x 0.008 = 0.02 pb. (158) 

The total is 

Acr qq - ~ 0.06 pb. (159) 

Since the perturbative cross section at m = 175 GeV is a^ 1 = 4.87, the extra contribution 

amounts to an increase of about 1 to 1.5%. 

For the gg-channel displayed in Fig. 8b, r) ~ 0.05, and we find 

Aa gg ~ \ (2 - 1) x 0.05 = 0.025 pb. (160) 

Given that a v gg rt ~ 0.65 pb, this correction is a 4% effect. We estimate that the total 
non-perturbative correction cannot be more than 

Aa tl ~ 0.02a% rt . (161) 

Correspondingly, we can expect a maximum cross section of about 5.7 pb at m = 175 GeV. 

7 Discussion and conclusions 

There are two main aspects of this paper, one of a more theoretical nature and the other 
phenomenological. First, we describe in all generality features of perturbative resummation 
that occur in a variety of hard scattering processes while focusing on the specific process of ti 
production. This process has the advantage of being physically very interesting in pQCD in 
that it probes our understanding of the theory when multiple physical scales are involved. It 
is characterized by an unquestionably large physical scale, m, and a partonic variable r\ that, 
in turn, is the ratio of two large scales. It is an ideal process in which one may examine the 
systematics of partonic interactions in more detail than previously possible. In this process, 
we encounter difficulties that have challenged the applicability of pQCD for a long time, 
notably IR renormalons and large color factors, as well as the systematics of universality 
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and invariance under changes of the hard scale. In dealing with tt production, we treat all 
these issues in a concrete, phenomenologically significant way. 

Second, the paper is concerned with predictions for the total tt cross section and the 
practical application of the theoretical resolutions developed in the first part. At the same 
time, comparisons are made with next-to-leading order cross sections and with earlier re- 
summation calculations, and extensive discussions are presented of theoretical limitations 
and uncertainties. 

Our theoretical analysis shows that perturbative resummation without a model for non- 
perturbative behavior is both possible and advantageous. In perturbative resummation, the 
perturbative region of phase space is separated cleanly from the region of non-perturbative 
behavior. The former is the region where large threshold corrections exponentiate but behave 
in a way that is perturbatively stable. The asymptotic behavior of the QCD perturbative 
series, including large multiplicative color factors, is flat, and excursions around the optimum 
number of perturbative terms does not create numerical instabilities or intolerable RG- 
dependence. Infrared renormalons are far away from the stability plateau and, even though 
their presence is essential for defining this plateau, they are of no numerical consequence 
in the perturbative regime. Large color factors, which are multiplicative, enhance the IR 
renormalon effects and contribute significantly to limiting the perturbative regime. 

As expected in a process with two physical scales, the constraint that insures good per- 
turbative behavior is a function of several parameters. Making the customary assumption 
that the renormalization scale is identified with the factorization scale (we denote this single 
scale by //), and denoting the variable that probes the partonic threshold by 77, we may 
summarize the essence of our perturbative resummation procedure as follows: 

The threshold logarithmic corrections are exponentiated in an exponent 



where x = ln(l/(l — z)). If the total center-of-mass energy \J~S is not far above threshold, 
i] ~ S/4m 2 — 1 ~ (1 — z)/2. The number of perturbative terms N(fj,) depends on the hard 
scale only, as long as IR renormalons are far away. The cross section is insensitive to fine- 
tuning of this number, as long as we remain in a specific well-defined perturbative region 
of x. In this paper, we determine this region in momentum space by demanding that the 
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cross section conform to perturbative power-counting. The cross section is a series multi- 
plying the exponentiated effects of the threshold and whose successive terms are given by 
functions of the successive derivatives d k E / k\dx k . Perturbative power-counting means that 
higher derivatives are suppressed relative to lower ones, because they contribute subleading 
logarithmic structures in ln(l/(l — z)). This statement of perturbative power-counting is 
the essence of perturbation theory. Imposing its validity in the two-scale parameter space, 
we determine the perturbative regime: 



This constraint restricts several sources of non-perturbative behavior, namely: IR renor- 
malons, x, and color factors. Some of them are channel- independent while others 

depend on the nature of the interacting partons. The constraint serves to identify the only 
unknown in the problem, i.e., the momentum range in x in which perturbative resummation 
is justified. 

We calculate 



Our total tt-production cross section at m = 175 GeV and y/S = 1.8 TeV is 10 to 14% 
greater than that of earlier calculations. Part of the increase comes from the more recent 
parton distributions we use. Our resummed cross sections are about 9% above the next-to- 
leading order cross sections computed with the same parton distributions. The renormaliza- 
tion/factorization scale dependence of our cross section is fairly flat, resulting in a 9 to 10% 
theoretical uncertainty. This variation is smaller than the corresponding dependence of the 
next-to-leading order cross section, and it is much smaller than the corresponding depen- 
dence of the resummed cross section of Ref . || . There are other perturbative uncertainties, 
such as dependence on parton distributions and factorization scheme. Each of these sources 
affects our cross section minimally, at level of 4% or less. These variations are not altogether 
independent, so we opt not to add them in estimating the theoretical uncertainty. For ex- 
ample, different parton distributions are associated with different values of A and therefore 
of a s . This uncertainty is correlated with uncertainty in a s from other measurements, and 
with the standard /i-variation that we use. Additional back-of-the-envelope modelling of 




(163) 



a u (m = 175 GeV, y/s = 1.8 TeV) = 5.52+°;°^ pb . 
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non-perturbative behavior does not increase our cross section more than 2%. However, we 
do not claim to address the issue of possible non-perturbative enhancements other than at 
the level of a conservative educated guess. 

Our theoretical analysis and the stability of our cross sections under \x variation provide 
confidence that our perturbative resummation procedure yields an accurate calculation of the 
inclusive top quark cross section at Tevatron energies and exhausts present understanding 
of the perturbative content of the theory. Our prediction agrees with data, within the large 
experimental uncertainties. We anticipate the greater precision of top quark production data 
from future Tevatron collider runs will be instrumental in a fundamental test of perturbative 
QCD. The data may also guide the modelling of the non-perturbative phase of the theory 
and provide a glimpse of phenomena beyond the Standard Model. 

The methodology described in this paper can be applied in several closely related sit- 
uations. The production of bottom quarks in hadron reactions at energies typical of the 
Fermilab fixed-target program or the HERA-B facility should be sensitive to the same type 
of threshold enhancements. Other reactions include the production of hadronic jets with 
very large values of transverse momentum or the production of very massive lepton pairs in 
Tevatron collider experiments. We hope to address these topics in the near future. 
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FIGURE CAPTIONS 

Figure 1. Optimum number of perturbative terms in the exponent. 

(a) Normalized principal-value exponents (solid) and their perturbative approxima- 
tions (dashed) as a function of N, for fixed mass and for four parametric mo- 
ment values (from the bottom, at n = 10,20,30,40). Optimization occurs at 
N(t(m = 175Gev)) = 6 for all four moments. 

(b) The function N(t(m)) (solid) and its simple analytic approximation [t(m) — 3/2] 
(dashed). 

Figure 2. The normalized principal-value exponent (solid) and its perturbative approxima- 
tion (dashed) versus the moment, n. 

Figure 3. The normalized perturbative exponent without regulators: 

(a) Partial sums versus the number N of terms as a parametric family in n: Solid 
bunch is for n G {100,1000} (increasing) in steps of 100. Dotted bunch is for 
n = 2000, 3000 (increasing). Dashed bunch is for n = 10000, 20000, 30000 
(increasing) . 

(b) Slope function versus N for n = 10, 50, 100, 500, 1000, 2000 (increasing). The 
minimum is attained at N — 1 ~ N(t(m = 175GeV)) = 6, for n < 1000 ~ m/A. 

Figure 4. Renormalization/factorization hard scale dependence of the resummation expo- 
nents: 
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(a) Drell-Yan perturbative (solid), principal value (dotted), and leading- logarithm 
(dashed) exponents versus /i for fixed mass and parametric moment values n = 
10, 50, 100 (increasing); N = N(t(m = 175GeV)) + 1. 

(b) Parametric families of partial sums N — 1, ...N(t(m = 175GeV)) + 1 (increasing) 
versus fi for the Drell-Yan (solid) and leading logarithm (dashed) exponents, at 
n = 50. 

Figure 5. Evaluation of the perturbative regime: 

(a) Saturation of the perturbative constraint in the qq channel as derived in moment 
space (S(n) = solid, a s e~ E ^ = dashed (Drell-Yan), dotted (leading logarithms)) 
and in momentum space (-Pi(n) < 1, dashed = Drell-Yan, dotted = leading loga- 
rithms); N = N(t(m = 175GeV)) + 1. The upper curves are those in momentum 
space. 

(b) Same as in (a) but for the gg channel. 

Figure 6. Hard-scale dependence of the perturbative regime: 

(a) The function P\(n) vs. n as a parametric family for /i = 100,150,200,250,300 
GeV for the Drell-Yan exponent (solid bunch) and the leading-logarithm exponent 
(dashed bunch, decreasing) for the qq channel. 

(b) Same as in (a) for the gg channel. 

Figure 7. The partonic cross sections a(rj) as a function of rj: 

(a) Born (dotted), next-to- leading order (dashed), and resummed (solid) partonic 
cross sections for the qq channel. 

(b) Same as in (a) for the gg channel. 

Figure 8. Differential cross section da/drj as a function of t] in the MS factorization scheme. 
The physical cross sections are the areas under the curves. 

(a) qq channel: Born (dotted), next-to-leading order (dashed), and resummed (solid). 
The non-perturbative regime is the area from 77 = to the point in r\ at which 
the solid and dashed curves intersect. 



(b) Same as in (a) for the gg channel. 



Figure 9. Inclusive cross section for top quark production in the MS scheme: 

(a) At the Tevatron, for pp — > ttX at v^S* =1.8 TeV. The extremum dashed lines are 
our band of uncertainty, and the solid line between them is our central- value pre- 
diction. We also reproduce the published data of the CDF and DO collaborations. 

(b) Hard-scale dependence of the resummed (solid) and next-to-leading order (dashed) 
cross sections at =1.8 TeV for m = 175 GeV. 

(c) Same as in (a), but for the upgraded Tevatron, v^S* = 2.0 TeV. 

(d) Central values of the resummed cross section (solid) for pp — > tiX at the CERN 
LHC energies of 10 and 14 TeV and the corresponding next-to-leading order pre- 
dictions (dashed). 

Figure 10. Physical cross sections in the qq channel in the MS scheme. The solid lines 
denote the finite-order partial sums of the universal leading-logarithmic contributions 
from the explicit 0(a) and 0(a 2 ) calculations for the ti and Drell-Yan cross sections, 
respectively, integrated throughout phase space. Lower solid: er^; middle solid: (j( 0+1 ); 
upper solid: a^ 0+1+2 \ The dashed curve represents the exact next-to-leading order 
calculation for tt production, in excellent agreement with a^ 0+1 ^. The dotted curve is 
our resummed prediction. 

Figure 11. Physical cross sections for the qq channel in the MS scheme: Leading-logarithm 
resummed cross sections for /j,/m — 1 (solid) and fi/m = 2 (dotted) and the Drell-Yan 
resummed version for (j,/m — 1 (dashed). 
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